*** COVID-19 economic and health consequence frames affect vaccine incentive attitudes in Germany – a framing experiment
*** Sebastian Jäckle & James K. Timmis

*** Version 16.05.25

clear all

*** install packages ***

	*ssc install tabcount, replace
	*ssc install mplotoffset, replace
	*ssc install blindschemes, replace
	*ssc install estout, replace
	*ssc install alphawgt, replace
	*ssc install catplot, replace
	*ssc install combomarginsplot, replace
	*ssc install coefplot, replace
	*net install grc1leg, from( http://www.stata.com/users/vwiggins/)
	*net install grc1leg2, from( http://digital.cgdev.org/doc/stata/MO/Misc)

set scheme plotplainblind

** set working directory
 cd "YOURPATH\replication"


******************************************************************
**																**
**		1)	Data Generation										**
**																**
******************************************************************

* Read data 
	use "data\PPD_July_2022_reduced.dta", replace


* Calculate Raking Weights based on gender, agegroup, state and highest education
	 egen age_group_raking = cut(v_125), at(18 , 31, 46, 61,  110)	 label 
	 table age_group_raking
	 label define age_group_raking 0 "18-30" 1 "31-45" 2 "46-60" 3 "> 60", modify
	 label values age_group_raking age_group_raking

* non-binary are coded as female for the raking
	 gen male_2 = v_370
	 replace male_2 = 2 if v_370 == 3 
	
* highest school leaving certificate
 	recode v_20 (1 2 9=1) (3 =2) (4 5 10 =3), gen(edu_3)
 	label define edu_3 1 "low" 2 "middle" 3 "high", modify
 	label values edu_3 edu_3
 	label variable edu_3					"educational attainment (low/middle/high)"
	
* Recode Sunday question for raking (CDU and CSU together)
	 recode v_1143 (2 8 = 1) (3 = 2) (4 = 3) (5 = 4) (6 = 5) (7 = 6)  (10/99 = 7), gen(sonntagsfrage)
	 label define sonntagsfrage 1 "CDU/CSU" 2 "SPD" 3 "Grüne" 4 "FDP" 5 "Linke" 6 "AfD" 7 "Sonstige", modify
	 label value sonntagsfrage sonntagsfrage

	 gen bundesland = v_19
	 replace bundesland =. if v_19 == 0
	 tab bundesland

* generate base weight = 1  
	gen wt = 1

 	* Vorgehen nach : https://www.stata-press.com/books/sw-preview.pdf ab S. 58
 	* wichtig! nicht die exakten Populationsgrößen verwenden, sondern gerundet auf tausend (dabei auf gleiche Anzahl an totals achten!)
 	* Mit zu großen Zahlen kommt svycal nicht zurecht


	* official data
 	#delimit ; 
 	svycal rake  i.bundesland i.age_group_raking  i.male_2 i.sonntagsfrage [pw=wt] , gen(rake_wt_2) ll(0.1) ul(10) tolerance(0.01) 
 	totals( _cons=61732
 			1.bundesland  =  7831
 			2.bundesland  =  9602
 			3.bundesland  =  2482
 			4.bundesland  =  2070
 			5.bundesland  =  472
 			6.bundesland  =  1293
 			7.bundesland  =  4422
 			8.bundesland  =  1332
 			9.bundesland  =  6145
 			10.bundesland  =  13104
 			11.bundesland  =  3091
 			12.bundesland  =  760
 			13.bundesland  =  3298
 			14.bundesland  =  1807
 			15.bundesland  =  2279
 			16.bundesland  =  1742
 			0.age_group_raking  =  9782
 			1.age_group_raking  =  12632
 			2.age_group_raking  =  16556
 			3.age_group_raking  =  22761
 			1.male_2  =  29870
 			2.male_2  =  31862
 			1.sonntagsfrage  =  16348
 			2.sonntagsfrage  =  12300
 			3.sonntagsfrage  =  14946
 			4.sonntagsfrage  =  4359
 			5.sonntagsfrage  =  2958
 			6.sonntagsfrage  =  6228
 			7.sonntagsfrage  =  4593
 			)	;
 #delimit cr	

 rename rake_wt_2 gewicht3

 label variable gewicht3 "Raking with state, age group, gender and voting intention"


 
*** Label variables

	* Sunday Question / Voting intention 
	label define sonntagsfrage 1 "CDU/CSU" 2 "SPD" 3 "Greens" 4 "FDP" 5 "Left-Party" 6 "AfD" 7 "Others", modify
	label value sonntagsfrage sonntagsfrage
	label variable sonntagsfrage "Voting intention"

	* Education
	label variable edu_3		"Educational attainment"

	* Solidarity
	gen solidarity_1 = v_1295
	label variable solidarity_1 "Those who pay a lot of taxes are already doing enough for society."
	gen solidarity_2 = v_1296
	label variable solidarity_2 "Society depends on the voluntary commitment of the population."
	gen solidarity_3 = v_1297
	label variable solidarity_3 "It is more important to work for the good of the community than for your own good."
	gen solidarity_4 = v_1298
	label variable solidarity_4 "If a stranger is in need of my support, I try to give it to them."
	gen solidarity_5 = v_1300
	label variable solidarity_5 "Important matters take precedence, even when free time beckons."
	label define solidarity 1 "completely disagree" 2 "rather disagree" 3 "neutral" 4 "rather agree" 5 "completely agree"
	label values solidarity_* solidarity

	* test solidarity index
	pwcorr solidarity_*
	factor solidarity_*
		rotate
		estat kmo
		alphawgt solidarity_* [aweight=gewicht3], item

	* --> solidarity index including solidarity_1 - solidarity_4 but not solidarity_5

	* reverse solidarity_1
	gen solidarity_1_rev = 6-solidarity_1

	* calculate solidarity index
	gen solidarity_index = (((solidarity_1_rev + solidarity_2 + solidarity_3 + solidarity_4)/4)-1)/4
	label variable solidarity_index "Solidarity"

	* Attitudes towards different school systems
	gen school_system_1 = (v_1301 - 1)/4
	label variable school_system_1 "Public/State schools"
	gen school_system_2 = (v_1302 - 1)/4
	label variable school_system_2 "Private schools"
	gen school_system_3 = (v_1303 - 1)/4
	label variable school_system_3 "Waldorf education"
	label define school_system 0 "0"  1 "1"
	label values school_system_* school_system

	* Attitudes towards different medical systems
	gen medical_system_1 = (v_1415 - 1)/4
	label variable medical_system_1 "Mainstream medicine"
	gen medical_system_2 = (v_1416 - 1)/4
	label variable medical_system_2 "Homeopathy"
	gen medical_system_3 = (v_1417 - 1)/4
	label variable medical_system_3 "Physiotherapy"
	gen medical_system_4 = (v_1418 - 1)/4
	label variable medical_system_4 "Acupuncture and acupressure"
	label define medical_system 0 "0"  1 "1"
	label values medical_system_* medical_system

	* left-right 0-1
	gen left_right = (v_773 -1)/10
	label define left_right 0 "left" 1 "right"
	label value left_right left_right
	label variable left_right "Left = 0 / Right = 1"

	* GAL-TAN
	gen GAL_TAN = (v_774 -1)/10
	label define GAL_TAN 0 "GAL" 1 "TAN"
	label value GAL_TAN GAL_TAN
	label variable GAL_TAN "GAL = 0 / TAN = 1"

	* Big-5
				* Big Five Inventory (BFI-10)
				*Autor/in: Rammstedt, B., Kemper, C. J., Klein, M. C., Beierlein, C., & Kovaleva, A.
				* In ZIS seit: 2014
				* DOI: https://doi.org/10.6102/zis76
				
				* Auswertungshinweise: Um Messwerte für die individuelle Ausprägung der Befragungsperson auf den fünf Persönlichkeitsdimensionen zu erhalten, 
				* werden pro Dimension die Antworten auf den beiden Items gemittelt. Hierzu wird zunächst das jeweils negativ gepolte Item rekodiert (Items 1, 3, 4, 5 und 7) 
				* und anschließend pro Dimension der Mittelwert aus dem rekodierten und dem nicht rekodierten Item gebildet. 
				* Der Wertebereich der fünf Dimensionen liegt dann jeweils zwischen 1 und 5 (für Referenzwerte siehe Tabelle 8 - Tabelle 12 oder Rammstedt, 2007).
	gen big_5_1 = v_1280
	label variable big_5_1 "I see myself as someone who is reserved."
	gen big_5_2 = v_1281
	label variable big_5_2 "I see myself as someone who is generally trusting."
	gen big_5_3 = v_1282
	label variable big_5_3 "I see myself as someone who tends to be lazy."
	gen big_5_4 = v_1283
	label variable big_5_4 "I see myself as someone who is relaxed, handles stress well."
	gen big_5_5 = v_1284
	label variable big_5_5 "I see myself as someone who has few artistic interests."
	gen big_5_6 = v_1285
	label variable big_5_6 "I see myself as someone who is outgoing, sociable."
	gen big_5_7 = v_1286
	label variable big_5_7 "I see myself as someone who tends to find fault with others."
	gen big_5_8 = v_1287
	label variable big_5_8 "I see myself as someone who does a thorough job."
	gen big_5_9 = v_1288
	label variable big_5_9 "I see myself as someone who gets nervous easily."
	gen big_5_10 = v_1289
	label variable big_5_10 "I see myself as someone who has an active imagination."
	label define big_5 1 "does not apply at all" 5 "fully applies"
	label values big_5_* big_5

	* recode negative poled items
	gen big_5_1_recode = 6 - big_5_1
	gen big_5_3_recode = 6 - big_5_3
	gen big_5_4_recode = 6 - big_5_4
	gen big_5_5_recode = 6 - big_5_5
	gen big_5_7_recode = 6 - big_5_7

	* Means for Big-5 Factors
		* Extraversion
			gen big_5_extraversion = (((big_5_1_recode + big_5_6)/2)-1)/4
		* Compatibility (Verträglichkeit)
			gen big_5_compatibility = (((big_5_2 + big_5_7_recode)/2)-1)/4
		* Conscientiousness (Gewissenhaftigkeit)
			gen big_5_conscientiousness = (((big_5_3_recode + big_5_8)/2)-1)/4
		* Neuroticism
			gen big_5_neuroticism = (((big_5_4_recode + big_5_9)/2)-1)/4
		* Openness for experience
			gen big_5_openness = (((big_5_5_recode + big_5_10)/2)-1)/4
			
			label variable big_5_extraversion "Extraversion (Big 5)"
			label variable big_5_neuroticism "Neuroticism (Big 5)"
			label variable big_5_openness "Openness (Big 5)"
			label variable big_5_conscientiousness "Conscientiousness (Big 5)"
			label variable big_5_compatibility "Compatibility (Big 5)"
			
			
	* Covid-19 Vaccination status
	gen vacc_status = v_1133
	label variable vacc_status "SARS-CoV-2 vaccination status"
	label define vacc_status 1 "not vaccinated" 2 "one vaccination received" 3 "two vaccinations received" 4 "three vaccinations received" 5 "four vaccinations received"
	label value vacc_status vacc_status	

	gen vacc_04 = vacc_status - 1
	label variable vacc_04 "# SARS-CoV-2 vaccine doses"
	label define vacc_04 0 "0" 1 "1" 2 "2" 3 "3" 4 "4"
	label value vacc_04 vacc_04	
	
	* vacc dummy
	gen vacc_dummy = .
	replace vacc_dummy = 0 if vacc_status  < 3
	replace vacc_dummy = 1 if vacc_status >= 3
	
	label variable vacc_dummy "SARS-CoV-2 Vaccination Dummy"
	label define vacc_dummy 0 "0 or 1" 1 ">=2"
	label value vacc_dummy vacc_dummy
		
		* Reason for getting vaccinated
	gen vacc_reason_1 = (v_1335 -1) /4
	label variable vacc_reason_1 "Protecting self"
	gen vacc_reason_2 = (v_1336 - 1)/4
	label variable vacc_reason_2 "Protecting others"
	gen vacc_reason_3 = (v_1337 - 1)/4
	label variable vacc_reason_3 "Participation in public events"
	gen vacc_reason_4 = (v_1338 - 1)/4
	label variable vacc_reason_4 "Vocational mandates"
	gen vacc_reason_5 = (v_1339 - 1)/4
	label variable vacc_reason_5 "Medical advice/recommendation of the Permanent Vaccination Commission (STIKO)"
	gen vacc_reason_6 = (v_1340 - 1)/4
	label variable vacc_reason_6 "Peer pressure"
	label define vacc_reason 0 "not at all important"  1 "very important"
	label value vacc_reason_* vacc_reason

	* How useful do you consider routine pediatric immunization, e.g. against measles/mumps/rubella (MMR)?  
	gen child_vaccinations = (v_1372 - 1) /3
	label variable child_vaccinations "How useful do you consider routine pediatric immunization, e.g. against measles/mumps/rubella (MMR)?"
	label define child_vaccinations 0 "not at all meaningful"  1 "very meaningful"
	label values child_vaccinations* child_vaccinations 

	* Voting intention (Sonntagsfrage)
	gen voting_intention_all = v_1143
	gen voting_intention_short = sonntagsfrage
	label values voting_intention_short sonntagsfrage

	* gender including non-binary
		gen gender = v_370

		* gender (including non-binary)
		label define gender 1 "male" 2 "female" 3 "non-binary" 
		label values gender gender	
		label variable gender "Gender"


	* Household income
	recode v_1131 (1=5) (2=4) (3=3) (4=2) (5=1) (8 = 6), gen(hh_income)
	label define hh_income 1 "I find it very difficult to make ends meet" 2 "I find it somewhat difficult to make ends meet" 3 "I can make ends meet" 4 "I can live somewhat comfortably on the income" 5 "I can live very comfortably on the income" 6 "does not apply (no own household income)"
	label value hh_income hh_income
	label variable hh_income "Household income"

	* Religious denomination
	gen religion = v_1311
	label variable religion "Religious denomination"
	label define religion 1 "non-denominational" 2 "Roman-Catholic" 3 "Protestant" 4 "Evangelical Free Church" 5 "Orthodox"  6 "Jewish" 7 "Muslim" 8 "Other"
	label value religion religion

	* Bundesland
	label define bundesland 1 "Baden-Wuerttemberg" 2 "Bavaria" 3 "Berlin" 4 "Brandenburg" 5 "Bremen" 6 "Hamburg" 7 "Hesse" 8 "Mecklenburg-Vorpommern" 9 "Lower Saxony" 10 "North Rhine-Westphalia" 11 "Rhineland Palatinate" 12 "Saarland" 13 "Saxony" 14 "Saxony-Anhalt"  15 "Schleswig Holstein" 16 "Thuringia" 17 "Foreign country"
	label values bundesland bundesland

	* age
	gen age = v_125

	* age group raking
	replace age_group_r = age_group_r + 1
	label define age_group_r 1 "18-30" 2 "31-45" 3 "46-60" 4 "> 60"
	label value age_group_r age_group_r
	label variable age_group_r "Age group"

	* weight(gewicht3) for youth = 1 and for bundesland = foreign = 1
	replace gewicht3 = 1 if age < 18
	replace gewicht3 = 1 if bundesland == 17

		* Covid-19 Threat
	gen threat_covid = (v_1374-1)/4
	label variable threat_covid "How threatening? COVID-19"
	label define threat_covid 0 "not at all threatening" 1 "very much threatening"
	label values threat_covid threat_covid

	* Spaltungsindex --> Korrigieren kenne nicht Kategorie auf missing setzen
			* don't know category --> missing
			gen v_1384a = v_1384
			replace v_1384a = . if v_1384 == 6

			gen v_1385a = v_1385
			replace v_1385a = . if v_1385 == 6

			gen v_1386a = v_1386
			replace v_1386a = . if v_1386 == 6

			gen v_1387a = v_1387
			replace v_1387a = . if v_1387 == 6

			gen v_1388a = v_1388
			replace v_1388a = . if v_1388 == 6

			gen v_1393a = v_1393
			replace v_1393a = . if v_1393 == 6
			
	factor v_1384a v_1385a v_1386a v_1387a v_1388a v_1393a
		rotate
		estat kmo
		alpha v_1384a v_1385a v_1386a v_1387a v_1388a v_1393a, item

	* --> besser ohne Ost/West und ohne Einkommenspaltung aber cancel-culture hat leider sehr viele missings

	factor v_1384a v_1385a v_1386a v_1387a v_1388a
		rotate
		estat kmo
		alpha v_1384a v_1385a v_1386a v_1387a v_1388a, item	
		
	gen spaltungsindex = (((v_1384a + v_1385a + v_1386a + v_1387a + v_1388a)/5)-1)/4
	label variable spaltungsindex "Perception of social division"

		* Decision to get vaccinated or to get not vaccinated...
	gen vacc_decision_1  = (v_1400-1)/4
	label variable vacc_decision_1 "... was very difficult and emotionally highly stressful for me."
	gen vacc_decision_2 = (v_1401-1)/4
	label variable vacc_decision_2 "... didn't take long for me."
	label define vacc_decision 0 "completely disagree" 1 "completely agree"
	label values vacc_decision_* vacc_decision

		* Current Covid-19 safety perception
	gen safety_perception = (-(v_1342-4))/3
	label variable safety_perception "Current COVID-19 safety perception"
	label define safety_perception 0 "very unsafe" 1 "very safe"
	label value safety_perception safety_perception

		* Faked Vaccination Certificate
	recode v_1346 (2 = 0) (1 = 1), gen(faked_certificate)
	label variable faked_certificate "Do you know anyone in your immediate circle who is using a fake vaccination certificate?"
	label define yes_no 0 "no" 1 "yes"
	label values faked_certificate yes_no

		* Talking about Covid-19 Vaccination
	recode v_1359 (2 = 0) (1 = 1), gen(conversation_vaccination)
	label variable conversation_vaccination "Do you occasionally have conversations with others about Corona vaccination?"
	label values conversation_vaccination yes_no	

		* Informed about Covid-19 Vaccination compared to conversation partners
		* "How well informed do you consider yourself to be on the subject of Corona vaccination compared to your interlocutors?"
	gen informed_vaccination = (v_1402-1)/4
	label variable informed_vaccination "Self perception knowledge COVID-19 vaccination"
	label define informed_vaccination 0 "much worse informed" 1 "much better informed"
	label values informed_vaccination informed_vaccination

		* With which of these people/groups are you reluctant to talk about Corona vaccination? 
	gen reluctant_family = v_1361
	gen reluctant_friends = v_1362
	gen reluctant_colleagues = v_1363
	gen reluctant_gp = v_1364
	gen reluctant_foreigners = v_1365
	gen reluctant_acquaintances = v_1433
	gen reluctant_noproblem = v_1426

		* Which Covid-19 protection measures meaningful
	gen protection_measures_1 = v_1348
	label variable protection_measures_1 "Distance regulations"
	gen protection_measures_2 = v_1349
	label variable protection_measures_2 "Access to certain areas only with negative test"
	gen protection_measures_3 = v_1350
	label variable protection_measures_3 "Wear protective mask in closed rooms"
	gen protection_measures_4 = v_1406
	label variable protection_measures_4 "none, all measures should be abolished."

		* Framing texts
	gen framing_long_covid_1 = v_1438
	gen framing_long_covid_2 = dupl1_v_1438
	gen framing_economic_1 = dupl2_v_1438
	gen framing_economic_2 = dupl3_v_1438

	gen framing = .
	replace framing = 0 if dispcode != 22 
	replace framing = 1 if framing_long_covid_1 == 1
	replace framing = 2 if framing_long_covid_2 == 1
	replace framing = 3 if framing_economic_1 == 1
	replace framing = 4 if framing_economic_2 == 1
	label define framing 0 "Control group" 1 "Long COVID frame 1" 2 "Long COVID frame 2" 3 "Economic frame 1 " 4 "Economic frame 2" 
	label value framing framing
	
	gen framing_pool_types = framing
	replace framing_pool_types = 1 if framing == 2
	replace framing_pool_types = 2 if framing == 3
	replace framing_pool_types = 2 if framing == 4
	label define framing_pool_types 0 "Control group" 1 "Long COVID frame" 2 "Economic frame"
	label value framing_pool_types framing_pool_types

	gen framing_pool_all = framing
	replace framing_pool_all = 1 if framing == 2
	replace framing_pool_all = 1 if framing == 3
	replace framing_pool_all = 1 if framing == 4
	label define framing_pool_all 0 "Control group" 1 "Framing group"
	label value framing_pool_all framing_pool_all	
	
	
		* Incentives for vaccination in general
	gen incentive_football = (v_1352-1)/4
	label variable incentive_football "Football ticket"
	gen incentive_sausage = (v_1353-1)/4
	label variable incentive_sausage "Free grilled sausage"
	gen incentive_lottery = (v_1354-1)/4
	label variable incentive_lottery  "Lottery"
	gen incentive_money = (v_1355-1)/4
	label variable incentive_money "Money"
	label define incentive 0 "not at all meaningful"  1 "very meaningful"
	label values incentive_* incentive

		* How much money as incentive for vaccination? Original ordinal scale
	gen money_incentive = v_1357
	label variable money_incentive "Provided there was a direct monetary compensation for Corona vaccination, how much money would you consider appropriate?"
	label define money_incentive 1 "nothing" 2 "up to 19.99 EUR" 3 "20 to 49,99 EUR" 4 "50 to 199.99 EUR" 5 "200 to 499.99 EUR" 6 "500 EUR and more"
	label values money_incentive money_incentive


		* How much money as incentive for vaccination? using mean values of categories
	gen money_incentive_metric = .
	replace money_incentive_metric = 0 		if money_incentive == 1
	replace money_incentive_metric = 10 	if money_incentive == 2
	replace money_incentive_metric = 35 	if money_incentive == 3
	replace money_incentive_metric = 125 	if money_incentive == 4
	replace money_incentive_metric = 350 	if money_incentive == 5
	replace money_incentive_metric = 500 	if money_incentive == 6

		* Participation in treatment costs for unvaccinated persons
	gen cost_share_yes = v_1395
	label variable cost_share_yes "If unvaccinated people fall ill with COVID-19 and have to go to hospital, should they share in the cost of treatment?"
	label define cost_share_yes 1 "yes" 2 "no" 3 "no opinion"
	label value cost_share_yes cost_share_yes
	gen cost_share_percentage = v_1396
	label variable cost_share_percentage "What percentage of treatment costs should unvaccinated people pay?"

		* filter for cost_share (yes or no, but not "no opinion"
	gen filter_cost_share = 0
	replace filter_cost_share = 1 if cost_share_yes < 3


		* Covid-19 compulsory vaccination
	gen vaccination_compulsory = v_1414
	label variable vaccination_compulsory "Are you in favour of compulsory Corona vaccination for this group of persons?"
	label define vaccination_compulsory 1 "for all adults" 2 "for all elderly and at-risk groups" 3 "for certain professions" 4 "no, I reject compulsory vaccination"

		* Compulsory vaccination against measles
	gen measles_compulsory = v_1371
	label variable measles_compulsory "Since March 2020, vaccination against measles has been compulsory in Germany for certain groups of people (children, teachers, care staff in kindergartens, etc.). What do you think about this?"
	label define measles 1 "good" 2 "bad" 3 "no opinion"


		* healthstate
	recode v_498 (1=5) (2=4) (3=3) (4=2) (5=1), gen(healthstate)
	label define healthstate 1 "very bad" 2 "bad" 3 "average" 4 "good"  5 "very good"
	label value healthstate healthstate
	label variable healthstate "Health status"

	gen health_01 = (healthstate-1)/4
	label define health_01 0 "very bad"   1 "very good"
	label value health_01 health_01
	label variable health_01 "Health status"

		* Migration background (yes/no)
	gen migration = v_57
	replace migration = 0 if v_57 == 2

		* Highest general school leaving certificate
	gen edu_1 = v_20
	label variable edu_1 "Highest general school leaving certificate"
	label define edu_1 1 "none" 2 "Volks-/ Hauptschulabschluss bzw. Polytechnische Oberschule mit Abschluss 8. oder 9. Klasse" 3 "Mittlere Reife, Realschulabschluss bzw. Polytechnische Oberschule mit Abschluss 10. Klasse" 4 "Fachabitur" 5 "Abitur bzw. Erweiterte Oberschule mit Abschluss 12. Klasse (Hochschulreife)" 9 "Other" 10 "still in school" 

		* Highest vocational or academic certificate
	gen edu_2 = v_1436
	label variable edu_2 "Highest vocational or academic certificate"
	label define edu_2 1 "none" 2 "completed apprenticeship" 3 "Technical school or vocational school diploma" 4 "Master craftsman/technician or equivalent technical college degree" 5 "(Technical) college degree (also degree from a school of engineering)" 6 "University degree" 7 "PhD" 8 "Other" 

		* Employment status
	gen fulltime = v_1123
	label variable fulltime "I work fulltime"
	gen parttime = v_1124
	label variable parttime "I work parttime"
	gen minijob = v_1125
	label variable minijob "Minijob"
	gen school = v_1126
	label variable school "I go to school"
	gen vocational_training = v_1127
	label variable vocational_training "I am doing vocational training"
	gen university = v_1128
	label variable university "I study at a university"
	gen parental_leave = v_1142
	label variable parental_leave "I am on parental leave or maternity leave at the moment"
	gen unemployed = v_1129
	label variable unemployed "I am unemployed"
	gen pension = v_1130
	label variable pension "I receive a pension"



		* Private health insurance
	recode v_1408 (1 = 0) (2 = 1), gen(private_health_insurance)
	label variable private_health_insurance "Private health insurance"


		

		* gewicht for youth = 1 and for bundesland = foreign = 1
	replace gewicht = 1 if age < 18
	replace gewicht = 1 if bundesland == 17

	**** during the first hours of the survey an identifier for the frame was missing. These cases have to be deleted.
	drop if startzeitpunkt < (tc(30jun2022 16:04:22))

	*** only cases that had reached the frames
	drop if framing == .

	*** Generate filter variables for listwise deletion
		
	* Filter for listwise deletion in Analysis I
	gen filter_1 = 1
	replace filter_1 = 0 if incentive_football == .
	replace filter_1 = 0 if incentive_sausage == .
	replace filter_1 = 0 if incentive_lottery == .
	replace filter_1 = 0 if incentive_money == .
	replace filter_1 = 0 if vacc_04 == .
	replace filter_1 = 0 if age_group_r == .
	replace filter_1 = 0 if gender == .
	replace filter_1 = 0 if left_right == .
	replace filter_1 = 0 if GAL_TAN == .
	replace filter_1 = 0 if solidarity_index == .
	replace filter_1 = 0 if big_5_e == .
	replace filter_1 = 0 if big_5_n == .
	replace filter_1 = 0 if big_5_com == .
	replace filter_1 = 0 if big_5_con == .
	replace filter_1 = 0 if big_5_o == .
	replace filter_1 = 0 if sonntagsfrage == .
	replace filter_1 = 0 if medical_system_1 == .
	replace filter_1 = 0 if medical_system_2 == .
	replace filter_1 = 0 if school_system_3 == .
	replace filter_1 = 0 if religion == .
	replace filter_1 = 0 if bundesland == .
	replace filter_1 = 0 if edu_3 == .
	replace filter_1 = 0 if hh_inc == .
	replace filter_1 = 0 if health_01 == .
	replace filter_1 = 0 if spaltungsindex == .
	replace filter_1 = 0 if threat_covid == .
	replace filter_1 = 0 if safety_perception == .
	replace filter_1 = 0 if vacc_decision_1 == .
	replace filter_1 = 0 if vacc_decision_2 == .
	replace filter_1 = 0 if money_incentive == .
	replace filter_1 = 0 if money_incentive_metric == .		
	
	* Save data 
	save "dump\PPD_July_2022_labeled.dta", replace
	
	

	
******************************************************************
**																**
**		3)	Regression Models									**
**																**
******************************************************************		
	use "dump\PPD_July_2022_labeled.dta", replace
	
	** ANALYSIS I: DV = Attitudeds towards incentives
	
	* 3.0_all) bivariate regression based on all cases for which the framing variable and the dependent variable (attitudes towards incentives) are available (no listwise deletion)
	ologit incentive_football i.framing
	est store m0_ologit_football_all
	reg  incentive_football i.framing 
	est store m0_ols_football_all
		* for those with at least 1 dose vs those without
		reg  incentive_football i.framing if vacc_status != 1
		reg  incentive_football i.framing if vacc_status == 1
	
	ologit incentive_sausage i.framing 
	est store m0_ologit_sausage_all
	reg incentive_sausage i.framing 
	est store m0_ols_sausage_all
		* for those with at least 1 dose vs those without
		reg  incentive_sausage i.framing if vacc_status != 1
		reg  incentive_sausage i.framing if vacc_status == 1

	ologit incentive_lottery i.framing 
	est store m0_ologit_lottery_all
	reg incentive_lottery i.framing 
	est store m0_ols_lottery_all
		* for those with at least 1 dose vs those without
		reg  incentive_lottery i.framing if vacc_status != 1
		reg  incentive_lottery i.framing if vacc_status == 1

	ologit incentive_money i.framing 
	est store m0_ologit_money_all
	reg incentive_money i.framing  
	est store m0_ols_money_all
		* for those with at least 1 dose vs those without
		reg  incentive_money i.framing if vacc_status != 1
		reg  incentive_money i.framing if vacc_status == 1
	
	
	ologit money_incentive i.framing
	est store m0_ologit_howmuchmoney_all
	reg money_incentive_metric i.framing
	est store m0_ols_howmuchmoney_all
	
	* 3.0_lw) bivariate regression based on listwise deletion
	ologit incentive_football i.framing if filter_1 == 1
	est store m0_ologit_football_lw
	reg  incentive_football i.framing if filter_1 == 1
	est store m0_ols_football_lw


	ologit incentive_sausage i.framing if filter_1 == 1
	est store m0_ologit_sausage_lw
	reg incentive_sausage i.framing if filter_1 == 1
	est store m0_ols_sausage_lw

	ologit incentive_lottery i.framing if filter_1 == 1
	est store m0_ologit_lottery_lw
	reg incentive_lottery i.framing if filter_1 == 1
	est store m0_ols_lottery_lw

	ologit incentive_money i.framing if filter_1 == 1
	est store m0_ologit_money_lw
	reg incentive_money i.framing  if filter_1 == 1
	est store m0_ols_money_lw
	
	ologit money_incentive i.framing if filter_1 == 1
	est store m0_ologit_howmuchmoney_lw
	reg money_incentive_metric i.framing if filter_1 == 1
	est store m0_ols_howmuchmoney_lw
	
	
	* 3.1_ols) full regression model (OLS)
	reg incentive_football  													/// Analysis I DV 
			i.framing															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m1_ols_football
	margins, at(framing=(0(1)4)  ) saving(dump\football, replace)
	marginsplot, recast(scatter) title("Football ticket") ytitle("Meaningfulness") yline (0.2236774) xtitle("") xlabel(,angle(45)) ylabel(.15(.05)0.35)
	graph save "dump/fig_football_ols.gph", replace
	
	reg incentive_sausage 													/// Analysis I DV 
			i.framing															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m1_ols_sausage
	margins, at(framing=(0(1)4)  ) saving(dump\sausage, replace)
	marginsplot, recast(scatter) title("Sausage") ytitle("Meaningfulness") yline (0.2064297) xtitle("") xlabel(,angle(45)) ylabel(.15(.05)0.35)	
	graph save "dump/fig_sausage_ols.gph", replace	
	
	reg incentive_lottery  													/// Analysis I DV 
			i.framing															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m1_ols_lottery
	margins, at(framing=(0(1)4)  ) saving(dump\lottery, replace)
	marginsplot, recast(scatter) title("Lottery") ytitle("Meaningfulness") yline (0.2774382) xtitle("") xlabel(,angle(45)) ylabel(.15(.05)0.35)
	graph save "dump/fig_lottery_ols.gph", replace	
	
	reg incentive_money 													/// Analysis I DV 
			i.framing															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m1_ols_money
	margins, at(framing=(0(1)4)  ) saving(dump\money, replace)
	marginsplot, recast(scatter) title("Money") ytitle("Meaningfulness") yline (0.1907216) xtitle("") xlabel(,angle(45)) ylabel(.15(.05)0.35)
	graph save "dump/fig_money_ols.gph", replace	

	
	* 3.1_ols) full regression model (OLS) - framing_pool_types
	reg incentive_football  													/// Analysis I DV 
			i.framing_pool_types															/// framing_pooled_types
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m1_ols_football_pooled2
	margins, at(framing=(0(1)2)  ) saving(dump\football_pooled2, replace)
	marginsplot, recast(scatter) title("Football ticket") ytitle("Meaningfulness") yline (0.2236774) xtitle("") xlabel(,angle(45)) ylabel(.15(.05)0.35)
	graph save "dump/fig_football_ols_pooled2.gph", replace
	
	reg incentive_sausage 													/// Analysis I DV 
			i.framing_pool_types															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m1_ols_sausage_pooled2
	margins, at(framing=(0(1)2)  ) saving(dump\sausage_pooled2, replace)
	marginsplot, recast(scatter) title("Sausage") ytitle("Meaningfulness") yline (0.2064297) xtitle("") xlabel(,angle(45)) ylabel(.15(.05)0.35)	
	graph save "dump/fig_sausage_ols_pooled2.gph", replace	
	
	reg incentive_lottery  													/// Analysis I DV 
			i.framing_pool_types															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m1_ols_lottery_pooled2
	margins, at(framing=(0(1)2)  ) saving(dump\lottery_pooled2, replace)
	marginsplot, recast(scatter) title("Lottery") ytitle("Meaningfulness") yline (0.2774382) xtitle("") xlabel(,angle(45)) ylabel(.15(.05)0.35)
	graph save "dump/fig_lottery_ols_pooled2.gph", replace	
	
	reg incentive_money 													/// Analysis I DV 
			i.framing_pool_types															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m1_ols_money_pooled2
	margins, at(framing=(0(1)2)  ) saving(dump\money_pooled2, replace)
	marginsplot, recast(scatter) title("Money") ytitle("Meaningfulness") yline (0.1907216) xtitle("") xlabel(,angle(45)) ylabel(.15(.05)0.35)
	graph save "dump/fig_money_ols_pooled2.gph", replace	
	
	
	* 3.1_ols) full regression model (OLS) - framing_pool_all
	reg incentive_football  													/// Analysis I DV 
			i.framing_pool_all															/// framing_pooled_types
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m1_ols_football_pooled3
	margins, at(framing=(0(1)1)  ) saving(dump\football_pooled3, replace)
	marginsplot, recast(scatter) title("Football ticket") ytitle("Meaningfulness") yline (0.2236774) xtitle("") xlabel(,angle(45)) ylabel(.15(.05)0.35)
	graph save "dump/fig_football_ols_pooled3.gph", replace
	
	reg incentive_sausage 													/// Analysis I DV 
			i.framing_pool_all															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m1_ols_sausage_pooled3
	margins, at(framing=(0(1)1)  ) saving(dump\sausage_pooled3, replace)
	marginsplot, recast(scatter) title("Sausage") ytitle("Meaningfulness") yline (0.2064297) xtitle("") xlabel(,angle(45)) ylabel(.15(.05)0.35)	
	graph save "dump/fig_sausage_ols_pooled3.gph", replace	
	
	reg incentive_lottery  													/// Analysis I DV 
			i.framing_pool_all															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m1_ols_lottery_pooled3
	margins, at(framing=(0(1)1)  ) saving(dump\lottery_pooled3, replace)
	marginsplot, recast(scatter) title("Lottery") ytitle("Meaningfulness") yline (0.2774382) xtitle("") xlabel(,angle(45)) ylabel(.15(.05)0.35)
	graph save "dump/fig_lottery_ols_pooled3.gph", replace	
	
	reg incentive_money 													/// Analysis I DV 
			i.framing_pool_all															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m1_ols_money_pooled3
	margins, at(framing=(0(1)1)  ) saving(dump\money_pooled3, replace)
	marginsplot, recast(scatter) title("Money") ytitle("Meaningfulness") yline (0.1907216) xtitle("") xlabel(,angle(45)) ylabel(.15(.05)0.35)
	graph save "dump/fig_money_ols_pooled3.gph", replace	
	
	
	
	
	* 3.1_ologit) full regression model (Ologit)
	ologit incentive_football  													/// Analysis I DV 
			i.framing															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1
	est store m1_ologit_football

	ologit incentive_sausage 													/// Analysis I DV 
			i.framing															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1
	est store m1_ologit_sausage

	ologit incentive_lottery  													/// Analysis I DV 
			i.framing															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1
	est store m1_ologit_lottery

	ologit incentive_money 													/// Analysis I DV 
			i.framing															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1
	est store m1_ologit_money

	
	
	* 3.1a_ols) check with drivers (OLS)
	reg incentive_football  													/// Analysis I DV 
			i.framing															/// frames
			vacc_reason*														/// Extra-control: concrete drivers for vaccination decision
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 &  vacc_status != 1
	est store m1a_ols_football

	reg incentive_sausage 													/// Analysis I DV 
			i.framing															/// frames
			vacc_reason*														/// Extra-control: concrete drivers for vaccination decision
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 &  vacc_status != 1
	est store m1a_ols_sausage

	reg incentive_lottery  													/// Analysis I DV 
			i.framing															/// frames
			vacc_reason*														/// Extra-control: concrete drivers for vaccination decision
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 &  vacc_status != 1
	est store m1a_ols_lottery

	reg incentive_money 													/// Analysis I DV 
			i.framing															/// frames
			vacc_reason*														/// Extra-control: concrete drivers for vaccination decision
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 &  vacc_status != 1
	est store m1a_ols_money
	

	

	* 3.1b_ols) check with self perceived knowledge (OLS)
	reg incentive_football  													/// Analysis I DV 
			i.framing															/// frames
			informed_vaccination												/// Extra-control: self perceived knowledge COVID-19 vaccination
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m1b_ols_football

	reg incentive_sausage 													/// Analysis I DV 
			i.framing															/// frames
			informed_vaccination												/// Extra-control: self perceived knowledge COVID-19 vaccination
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m1b_ols_sausage

	reg incentive_lottery  													/// Analysis I DV 
			i.framing															/// frames
			informed_vaccination												/// Extra-control: self perceived knowledge COVID-19 vaccination
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m1b_ols_lottery

	reg incentive_money 													/// Analysis I DV 
			i.framing															/// frames
			informed_vaccination												/// Extra-control: self perceived knowledge COVID-19 vaccination
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m1b_ols_money
	
	
	

	* 3.1c_ols) check with private health insurance (OLS)
	reg incentive_football  													/// Analysis I DV 
			i.framing															/// frames
			private_health_insurance											/// Extra-control: private health insurance
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m1c_ols_football

	reg incentive_sausage 													/// Analysis I DV 
			i.framing															/// frames
			private_health_insurance											/// Extra-control: private health insurance
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m1c_ols_sausage

	reg incentive_lottery  													/// Analysis I DV 
			i.framing															/// frames
			private_health_insurance											/// Extra-control: private health insurance
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m1c_ols_lottery

	reg incentive_money 													/// Analysis I DV 
			i.framing															/// frames
			private_health_insurance											/// Extra-control: private health insurance
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m1c_ols_money
	
	
	
	* 3.2_ols)Interactions
	
	*** Interaction framing x age_group

	reg incentive_football  													/// Analysis I DV 
			i.framing##i.age_group_r											/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_football_int_age
	margins , at(age_group_r=(1(1)4) framing=(0(1)4)  ) 
	*marginsplot, title("Football ticket") legend(row(1)) xtitle("")
	mplotoffset, offset(0.1) recast(scatter) title("Football ticket") ytitle("")  legend(off) xtitle("")
	graph save "dump/fig_int_age_a.gph", replace
	
	reg incentive_sausage 													/// Analysis I DV 
			i.framing##i.age_group_r											/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_sausage_int_age
	margins , at(age_group_r=(1(1)4) framing=(0(1)4)  ) 
	*marginsplot, title("Sausage") xtitle("")
	mplotoffset, offset(0.1) recast(scatter) title("Sausage") ytitle("") legend(off) xtitle("")
	graph save "dump/fig_int_age_b.gph", replace
	
	reg incentive_lottery  													/// Analysis I DV 
			i.framing##i.age_group_r											/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_lottery_int_age
	margins , at(age_group_r=(1(1)4) framing=(0(1)4)  ) 
	*marginsplot, title("Lottery") xtitle("")
	mplotoffset, offset(0.1) recast(scatter) title("Lottery") ytitle("") legend(off) xtitle("")
	graph save "dump/fig_int_age_c.gph", replace
	
	reg incentive_money 													/// Analysis I DV 
			i.framing##i.age_group_r											/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_money_int_age
	margins , at(age_group_r=(1(1)4) framing=(0(1)4)  ) 
	*marginsplot, title("Money") xtitle("")
	mplotoffset, offset(0.1) recast(scatter) title("Money") ytitle("") legend(off) xtitle("")	
	graph save "dump/fig_int_age_d.gph", replace

	
		*** Interaction framing x gender

	reg incentive_football  													/// Analysis I DV 
			i.framing##i.gender													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_football_int_gender
	margins , at(gender=(1(1)2) framing=(0(1)4)  )
	mplotoffset, title("") ytitle("") legend(off) offset(0.1) recast(scatter) xtitle("")
	graph save "dump/fig_int_gender_a.gph", replace
	
	reg incentive_sausage 														/// Analysis I DV 
			i.framing##i.gender													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_sausage_int_gender
	margins , at(gender=(1(1)2) framing=(0(1)4)  )
	mplotoffset, title("") ytitle("") legend(off) offset(0.1) recast(scatter) xtitle("")
	graph save "dump/fig_int_gender_b.gph", replace
	
	reg incentive_lottery  														/// Analysis I DV 
			i.framing##i.gender													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_lottery_int_gender
	margins , at(gender=(1(1)2) framing=(0(1)4)  )
	mplotoffset, title("") ytitle("") legend(off) offset(0.1) recast(scatter) xtitle("")
	graph save "dump/fig_int_gender_c.gph", replace
	
	reg incentive_money 														/// Analysis I DV 
			i.framing##i.gender													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_money_int_gender
	margins , at(gender=(1(1)2) framing=(0(1)4)  )
	mplotoffset, title("") ytitle("") legend(off) offset(0.1) recast(scatter) xtitle("")
	graph save "dump/fig_int_gender_d.gph", replace
	
	
		*** Interaction framing x voting intention

	reg incentive_football  													/// Analysis I DV 
			i.framing##i.sonntagsfrage												/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_football_int_vote
	margins , at(sonntagsfrage=(1(1)7) framing=(0(1)4)  ) 
	mplotoffset, title("Football ticket") legend(row(1)) offset(0.1) recast(scatter) xtitle("")
	graph save "dump/fig_int_sonntagsfrage_a.gph", replace
	
	reg incentive_sausage 														/// Analysis I DV 
			i.framing##i.sonntagsfrage													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_sausage_int_vote
	margins , at(sonntagsfrage=(1(1)7) framing=(0(1)4)  ) 
	mplotoffset, title("Sausage") legend(row(1)) offset(0.1) recast(scatter) xtitle("")
	graph save "dump/fig_int_sonntagsfrage_b.gph", replace
	
	reg incentive_lottery  														/// Analysis I DV 
			i.framing##i.sonntagsfrage													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_lottery_int_vote
	margins , at(sonntagsfrage=(1(1)7) framing=(0(1)4)  ) 
	mplotoffset, title("Lottery") legend(row(1)) offset(0.1) recast(scatter) xtitle("")
	graph save "dump/fig_int_sonntagsfrage_c.gph", replace
	
	reg incentive_money 														/// Analysis I DV 
			i.framing##i.sonntagsfrage													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_money_int_vote
	margins , at(sonntagsfrage=(1(1)7) framing=(0(1)4)  ) 
	mplotoffset, title("Money") legend(row(1)) offset(0.1) recast(scatter) xtitle("")
	graph save "dump/fig_int_sonntagsfrage_d.gph", replace
	
	*** Interaction framing x solidarity index

	reg incentive_football  													/// Analysis I DV 
			i.framing##c.solidarity_index											/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_football_int_solid
	margins , at(solidarity_index=(0(0.1)1) framing=(0(1)4)  ) 
	marginsplot, title("Football ticket") legend(row(1)) xtitle("")
	graph save "dump/fig_int_solidarity_a.gph", replace
	
	reg incentive_sausage 													/// Analysis I DV 
			i.framing##c.solidarity_index											/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_sausage_int_solid
	margins , at(solidarity_index=(0(0.1)1) framing=(0(1)4)  ) 
	marginsplot, title("Sausage") xtitle("")
	graph save "dump/fig_int_solidarity_b.gph", replace
	
	reg incentive_lottery  													/// Analysis I DV 
			i.framing##c.solidarity_index											/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_lottery_int_solid
	margins , at(solidarity_index=(0(0.1)1) framing=(0(1)4)  ) 
	marginsplot, title("Lottery") xtitle("")
	graph save "dump/fig_int_solidarity_c.gph", replace
	
	reg incentive_money 													/// Analysis I DV 
			i.framing##c.solidarity_index											/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_money_int_solid
	margins , at(solidarity_index=(0(0.1)1) framing=(0(1)4)  ) 
	marginsplot, title("Money") xtitle("")
	graph save "dump/fig_int_solidarity_d.gph", replace	

		*** Interaction framing x vacc_status

	reg incentive_football  													/// Analysis I DV 
			i.framing##i.vacc_04													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_football_int_vacc_04
	margins , at(vacc_04=(0(1)4) framing=(0(1)4)  )
	mplotoffset, title("") ytitle("") legend(row(1))  xtitle("")  offset(0.1) recast(scatter) 
	graph save "dump/fig_int_vacc_04_a.gph", replace
	
	reg incentive_sausage 														/// Analysis I DV 
			i.framing##i.vacc_04													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_sausage_int_vacc_04
	margins , at(vacc_04=(0(1)4) framing=(0(1)4)  )
	mplotoffset, title("") ytitle("") legend(row(1))  xtitle("")  offset(0.1) recast(scatter) 
	graph save "dump/fig_int_vacc_04_b.gph", replace
	
	reg incentive_lottery  														/// Analysis I DV 
			i.framing##i.vacc_04													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_lottery_int_vacc_04
	margins , at(vacc_04=(0(1)4) framing=(0(1)4)  )
	mplotoffset, title("") ytitle("") legend(row(1))  xtitle("")  offset(0.1) recast(scatter)
	graph save "dump/fig_int_vacc_04_c.gph", replace
	
	reg incentive_money 														/// Analysis I DV 
			i.framing##i.vacc_04													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_money_int_vacc_04
	margins , at(vacc_04=(0(1)4) framing=(0(1)4)  )
	mplotoffset, title("") ytitle("") legend(row(1))  xtitle("")  offset(0.1) recast(scatter)
	graph save "dump/fig_int_vacc_04_d.gph", replace

		*** Interaction framing x health

	reg incentive_football  													/// Analysis I DV 
			i.framing##c.health_01													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_football_int_health
	margins , at(health_01=(0(0.25)1) framing=(0(1)4)  )
	marginsplot, title("Football ticket") legend(row(1))  xtitle("")
	graph save "dump/fig_int_health_a.gph", replace
	
	reg incentive_sausage 														/// Analysis I DV 
			i.framing##c.health_01													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_sausage_int_health
	margins , at(health_01=(0(0.25)1) framing=(0(1)4)  )
	marginsplot, title("Sausage") legend(row(1))  xtitle("")
	graph save "dump/fig_int_health_b.gph", replace
	
	reg incentive_lottery  														/// Analysis I DV 
			i.framing##c.health_01													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_lottery_int_health
	margins , at(health_01=(0(0.25)1) framing=(0(1)4)  )
	marginsplot, title("Lottery") legend(row(1))  xtitle("")
	graph save "dump/fig_int_health_c.gph", replace
	
	reg incentive_money 														/// Analysis I DV 
			i.framing##c.health_01													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_money_int_health
	margins , at(health_01=(0(0.25)1) framing=(0(1)4)  )
	marginsplot, title("Money") legend(row(1))  xtitle("")
	graph save "dump/fig_int_health_d.gph", replace

	
	*** Interaction framing x homeopathy

	reg incentive_football  													/// Analysis I DV 
			i.framing##c.medical_system_2													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_football_int_homeop
	margins , at(medical_system_2=(0(0.25)1) framing=(0(1)4)  )
	marginsplot, title("Football ticket") legend(row(1))  
	graph save "dump/fig_int_homeopathy_a.gph", replace
	
	reg incentive_sausage 														/// Analysis I DV 
			i.framing##c.medical_system_2													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_sausage_int_homeop
	margins , at(medical_system_2=(0(0.25)1) framing=(0(1)4)  )
	marginsplot, title("Sausage") legend(row(1))  
	graph save "dump/fig_int_homeopathy_b.gph", replace
	
	reg incentive_lottery  														/// Analysis I DV 
			i.framing##c.medical_system_2													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_lottery_int_homeop
	margins , at(medical_system_2=(0(0.25)1) framing=(0(1)4)  )
	marginsplot, title("Lottery") legend(row(1))  
	graph save "dump/fig_int_homeopathy_c.gph", replace
	
	reg incentive_money 														/// Analysis I DV 
			i.framing##c.medical_system_2													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_money_int_homeop
	margins , at(medical_system_2=(0(0.25)1) framing=(0(1)4)  )
	marginsplot, title("Money") legend(row(1))  
	graph save "dump/fig_int_homeopathy_d.gph", replace	
	
	*** Interaction framing x waldorf

	reg incentive_football  													/// Analysis I DV 
			i.framing##c.school_system_3													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_football_int_waldorf
	margins , at(school_system_3=(0(0.25)1) framing=(0(1)4)  )
	marginsplot, title("Football ticket") legend(row(1))  
	graph save "dump/fig_int_waldorf_a.gph", replace
	
	reg incentive_sausage 														/// Analysis I DV 
			i.framing##c.school_system_3													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_sausage_int_waldorf
	margins , at(school_system_3=(0(0.25)1) framing=(0(1)4)  )
	marginsplot, title("Sausage") legend(row(1))  
	graph save "dump/fig_int_waldorf_b.gph", replace
	
	reg incentive_lottery  														/// Analysis I DV 
			i.framing##c.school_system_3													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_lottery_int_waldorf
	margins , at(school_system_3=(0(0.25)1) framing=(0(1)4)  )
	marginsplot, title("Lottery") legend(row(1))  
	graph save "dump/fig_int_waldorf_c.gph", replace
	
	reg incentive_money 														/// Analysis I DV 
			i.framing##c.school_system_3													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_money_int_waldorf
	margins , at(school_system_3=(0(0.25)1) framing=(0(1)4)  )
	marginsplot, title("Money") legend(row(1))  
	graph save "dump/fig_int_waldorf_d.gph", replace	
	
	
	*** Interaction framing x mainstream medicine

	reg incentive_football  													/// Analysis I DV 
			i.framing##c.medical_system_1													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_football_int_mainmed
	margins , at(medical_system_1=(0(0.25)1) framing=(0(1)4)  )
	marginsplot, title("Football ticket") legend(row(1))  
	graph save "dump/fig_int_mainmed_a.gph", replace
	
	reg incentive_sausage 														/// Analysis I DV 
			i.framing##c.medical_system_1													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_sausage_int_mainmed
	margins , at(medical_system_1=(0(0.25)1) framing=(0(1)4)  )
	marginsplot, title("Sausage") legend(row(1))  
	graph save "dump/fig_int_mainmed_b.gph", replace
	
	reg incentive_lottery  														/// Analysis I DV 
			i.framing##c.medical_system_1													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_lottery_int_mainmed
	margins , at(medical_system_1=(0(0.25)1) framing=(0(1)4)  )
	marginsplot, title("Lottery") legend(row(1))  
	graph save "dump/fig_int_mainmed_c.gph", replace
	
	reg incentive_money 														/// Analysis I DV 
			i.framing##c.medical_system_1													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_money_int_mainmed
	margins , at(medical_system_1=(0(0.25)1) framing=(0(1)4)  )
	marginsplot, title("Money") legend(row(1))  
	graph save "dump/fig_int_mainmed_d.gph", replace		


*** Interaction framing x conscientiousness

	reg incentive_football  													/// Analysis I DV 
			i.framing##c.big_5_con											/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_football_int_big5con
	margins , at(big_5_con=(0(0.1)1) framing=(0(1)4)  ) 
	marginsplot, title("Football ticket") legend(row(1)) xtitle("")
	graph save "dump/fig_int_big_5_con_a.gph", replace
	
	reg incentive_sausage 													/// Analysis I DV 
			i.framing##c.big_5_con											/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_sausage_int_big5con
	margins , at(big_5_con=(0(0.1)1) framing=(0(1)4)  ) 
	marginsplot, title("Sausage") xtitle("")
	graph save "dump/fig_int_big_5_con_b.gph", replace
	
	reg incentive_lottery  													/// Analysis I DV 
			i.framing##c.big_5_con											/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_lottery_int_big5con
	margins , at(big_5_con=(0(0.1)1) framing=(0(1)4)  ) 
	marginsplot, title("Lottery") xtitle("")
	graph save "dump/fig_int_big_5_con_c.gph", replace
	
	reg incentive_money 													/// Analysis I DV 
			i.framing##c.big_5_con											/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_money_int_big5con
	margins , at(big_5_con=(0(0.1)1) framing=(0(1)4)  ) 
	marginsplot, title("Money") xtitle("")
	graph save "dump/fig_int_big_5_con_d.gph", replace	


*** Interaction framing x perception of social division

	reg incentive_football  													/// Analysis I DV 
			i.framing##c.spaltungsindex											/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_football_int_spalt
	margins , at(spaltungsindex=(0(0.1)1) framing=(0(1)4)  ) 
	marginsplot, title("Football ticket") legend(row(1)) xtitle("")
	graph save "dump/fig_int_spaltungsindex_a.gph", replace
	
	reg incentive_sausage 													/// Analysis I DV 
			i.framing##c.spaltungsindex											/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_sausage_int_spalt
	margins , at(spaltungsindex=(0(0.1)1) framing=(0(1)4)  ) 
	marginsplot, title("Sausage") xtitle("")
	graph save "dump/fig_int_spaltungsindex_b.gph", replace
	
	reg incentive_lottery  													/// Analysis I DV 
			i.framing##c.spaltungsindex											/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_lottery_int_spalt
	margins , at(spaltungsindex=(0(0.1)1) framing=(0(1)4)  ) 
	marginsplot, title("Lottery") xtitle("")
	graph save "dump/fig_int_spaltungsindex_c.gph", replace
	
	reg incentive_money 													/// Analysis I DV 
			i.framing##c.spaltungsindex											/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m3_ols_money_int_spalt
	margins , at(spaltungsindex=(0(0.1)1) framing=(0(1)4)  ) 
	marginsplot, title("Money") xtitle("")
	graph save "dump/fig_int_spaltungsindex_d.gph", replace	


** ANALYSIS II: DV = how much money for incentivation
	reg money_incentive_metric 													/// Analysis II DV 
			i.framing															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m2_ols_howmuchmoney	
	
	reg money_incentive_metric 													/// Analysis II DV 
			i.framing_pool_types															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m2_ols_howmuchmoney_pooled2	

	reg money_incentive_metric 													/// Analysis II DV 
			i.framing_pool_all															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m2_ols_howmuchmoney_pooled3		
	
	
	ologit money_incentive_metric 												/// Analysis II DV 
			i.framing															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m2_ologit_howmuchmoney	

	** check with drivers
	reg money_incentive_metric 													/// Analysis II DV 
			i.framing															/// frames
			vacc_reason*														/// Extra-control: concrete drivers for vaccination decision
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 &  vacc_status != 1
	est store m2a_ols_howmuchmoney	

	** check with self perceived knowledge
	reg money_incentive_metric 													/// Analysis II DV 
			i.framing															/// frames
			informed_vaccination												/// Extra-control: self perceived knowledge COVID-19 vaccination
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m2b_ols_howmuchmoney	

	** check with private health insurance
	reg money_incentive_metric 													/// Analysis II DV 
			i.framing															/// frames
			private_health_insurance											/// Extra-control: private health insurance
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m2c_ols_howmuchmoney	
	
*** Interactions
	*** Interaction framing x age_group

	reg money_incentive_metric 													/// Analysis II DV  
			i.framing##i.age_group_r											/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m4_ols_int_age
	margins , at(age_group_r=(1(1)4) framing=(0(1)4)  ) 
	*marginsplot, title("Age") legend(row(1)) xtitle("") ytitle("Money paid for vaccination") 
	mplotoffset, title("Age") legend(row(1)) offset(0.1) recast(scatter) xtitle("") ytitle("Money paid for vaccination")
	graph save "dump/fig_int2_age.gph", replace
	
	*** Interaction framing x gender

	reg money_incentive_metric 													/// Analysis II DV  
			i.framing##i.gender													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m4_ols_int_gender
	margins , at(gender=(1(1)2) framing=(0(1)4)  )
	mplotoffset, title("Gender") legend(row(1)) offset(0.1) recast(scatter) xtitle("") ytitle("Money paid for vaccination")
	graph save "dump/fig_int2_gender.gph", replace
	
	*** Interaction framing x voting intention

	reg money_incentive_metric 													/// Analysis II DV  
	i.framing##i.sonntagsfrage												/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m4_ols_int_vote
	margins , at(sonntagsfrage=(1(1)7) framing=(0(1)4)  ) 
	mplotoffset, title("Voting intention") legend(row(1)) offset(0.1) recast(scatter) xtitle("") ytitle("Money paid for vaccination")
	graph save "dump/fig_int2_sonntagsfrage.gph", replace	
	
	*** Interaction framing x solidarity index

	reg money_incentive_metric 													/// Analysis II DV  
			i.framing##c.solidarity_index											/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m4_ols_int_solid
	margins , at(solidarity_index=(0(0.1)1) framing=(0(1)4)  ) 
	marginsplot, title("Solidarity index") legend(row(1)) xtitle("") ytitle("Money paid for vaccination")
	graph save "dump/fig_int2_solidarity.gph", replace	
	
	*** Interaction framing x vacc_status

	reg money_incentive_metric 													/// Analysis II DV  
			i.framing##i.vacc_04													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m4_ols_int_vacc_04
	margins , at(vacc_04=(0(1)4) framing=(0(1)4)  )
	mplotoffset, title("How much money as incentive") legend(row(1))  xtitle("")  offset(0.1) recast(scatter) 
	graph save "dump/fig_int2_vacc_04.gph", replace	
	
	*** Interaction framing x health

	reg money_incentive_metric 													/// Analysis II DV  
			i.framing##c.health_01													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m4_ols_int_health
	margins , at(health_01=(0(0.25)1) framing=(0(1)4)  )
	marginsplot, title("How much money as incentive") legend(row(1))  xtitle("")
	graph save "dump/fig_int2_health.gph", replace

	*** Interaction framing x homeopathy

	reg money_incentive_metric 													/// Analysis II DV  
			i.framing##c.medical_system_2													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m4_ols_int_homeop
	margins , at(medical_system_2=(0(0.25)1) framing=(0(1)4)  )
	marginsplot, title("How much money as incentive") legend(row(1))  
	graph save "dump/fig_int2_homeopathy.gph", replace	

	*** Interaction framing x waldorf

	reg money_incentive_metric 													/// Analysis II DV  
			i.framing##c.school_system_3													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m4_ols_int_waldorf
	margins , at(school_system_3=(0(0.25)1) framing=(0(1)4)  )
	marginsplot, title("How much money as incentive") legend(row(1))  
	graph save "dump/fig_int2_waldorf.gph", replace	

	*** Interaction framing x mainstream medicine

	reg money_incentive_metric 													/// Analysis II DV  
			i.framing##c.medical_system_1													/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m4_ols_int_mainmed
	margins , at(medical_system_1=(0(0.25)1) framing=(0(1)4)  )
	marginsplot, title("How much money as incentive") legend(row(1))  
	graph save "dump/fig_int2_mainmed.gph", replace	

	*** Interaction framing x conscientiousness

	reg money_incentive_metric 													/// Analysis II DV  
			i.framing##c.big_5_con											/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m4_ols_int_big5con
	margins , at(big_5_con=(0(0.1)1) framing=(0(1)4)  ) 
	marginsplot, title("How much money as incentive") legend(row(1)) xtitle("")
	graph save "dump/fig_int2_big_5_con.gph", replace	
	
*** Interaction framing x perception of social division

	reg money_incentive_metric 													/// Analysis II DV  
			i.framing##c.spaltungsindex											/// Interaction with frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	est store m4_ols_int_spalt
	margins , at(spaltungsindex=(0(0.1)1) framing=(0(1)4)  ) 
	marginsplot, title("How much money as incentive") legend(row(1)) xtitle("")
	graph save "dump/fig_int2_spaltungsindex.gph", replace

	
	********************************************************************************************************************************************************************	
	* Table S2: Outcome I: Main models, OLS regression (b-coefficients)
	********************************************************************************************************************************************************************

	esttab m1_ols_football m1_ols_sausage m1_ols_lottery m1_ols_money  using "tables\table_S02.rtf", label replace r2 ar2 


	********************************************************************************************************************************************************************	
	* Table S3: Outcome I: Cross check models A) concrete drivers , OLS regression (b-coefficients)
	********************************************************************************************************************************************************************

	esttab m1a_ols_football m1a_ols_sausage m1a_ols_lottery m1a_ols_money  using "tables\table_S03.rtf", label replace r2 ar2 
	
	********************************************************************************************************************************************************************	
	* Table S4: Outcome I: Cross check models B) self perceived knowledge SARS-CoV-2 vaccination , OLS regression (b-coefficients)
	********************************************************************************************************************************************************************

	esttab m1b_ols_football m1b_ols_sausage m1b_ols_lottery m1b_ols_money  using "tables\table_S04.rtf", label replace r2 ar2 
	
	********************************************************************************************************************************************************************	
	* Table S5: Outcome I: Cross check models C) private health insurance , OLS regression (b-coefficients)
	********************************************************************************************************************************************************************

	esttab m1c_ols_football m1c_ols_sausage m1c_ols_lottery m1c_ols_money  using "tables\table_S05.rtf", label replace r2 ar2 
	
	********************************************************************************************************************************************************************	
	* Table S6: Outcome II: Main model + cross check models, OLS regression (b-coefficients)
	********************************************************************************************************************************************************************

	esttab m2_ols_howmuchmoney  m2a_ols_howmuchmoney	 m2b_ols_howmuchmoney	 m2c_ols_howmuchmoney	 using "tables\table_S06.rtf", label replace r2 ar2
	
******************************************************************
**																**
**		4)	Figures Main Article								**
**																**
******************************************************************	

	*****************************************************************
	* Fig 01: Participant flow chart		
	*****************************************************************	

	* numbers for flow chart
	
	* survey started
		use "data\PPD_July_2022.dta", replace
		**** during the first hours of the survey an identifier for the frame was missing. These cases have to be deleted.
		drop if startzeitpunkt < (tc(30jun2022 16:04:22))
		tab dispcode

	* Dropout before framing experiment & number in experiment and control groups
				* Framing texts
				gen framing_long_covid_1 = v_1438
				gen framing_long_covid_2 = dupl1_v_1438
				gen framing_economic_1 = dupl2_v_1438
				gen framing_economic_2 = dupl3_v_1438

				gen framing = .
				replace framing = 0 if dispcode != 22 
				replace framing = 1 if framing_long_covid_1 == 1
				replace framing = 2 if framing_long_covid_2 == 1
				replace framing = 3 if framing_economic_1 == 1
				replace framing = 4 if framing_economic_2 == 1
				label define framing 0 "Control group" 1 "Long COVID frame 1" 2 "Long COVID frame 2" 3 "Economic frame 1 " 4 "Economic frame 2" 
				label value framing framing
				
				tab framing, miss
		
	* framing experiment & item non response	
		drop if framing == .
		tab dispcode
	
	use "dump\PPD_July_2022_labeled.dta", replace
	
	* participants < 18
	tab age if age < 18
	
	* item non-response
	drop if age < 18
	tab filter_1
	
	
	
	*****************************************************************
	* Fig 02: Descriptive statistics of dependent variables by vaccination status		
	*****************************************************************
	
	use "dump\PPD_July_2022_labeled.dta", replace
	
	catplot incentive_football vacc_dummy if filter_1 == 1 , percent(vacc_dummy) stack asyvars recast(bar) ///
		fxsize(25)									///
		title("Free ticket for" "football match")	///
		ytitle("")									///
		l1title("", )								///
		ylabel(,)					///
		blabel(bar, pos(center) format(%2.1f) )	///
		var1opts(gap(*0) label(labsize(large))) ///
			bar(1, color(538t) fintensity(inten90)) ///
			bar(2, color(538t) fintensity(inten50)) ///
			bar(3, color(gray) fintensity(inten40)) ///
			bar(4, color(538o) fintensity(inten50)) ///
			bar(5, color(538o) fintensity(inten90)) ///
		legend(rows(1) stack size(vsmall) 	///
		symplacement(center) pos(6) 		///
		label(1 "not at all meaningful") 	///
		label(2 "rather not meaningful")	///
		label(3 "neutral") 	///
		label(4 "rather meaningful") 	///
		label(5 "very meaningful")) 	///
		text(50 -42 "percent", orientation(vertical) size(medsmall) )
		graph save "dump/fig_02a.gph", replace

	catplot incentive_sausage vacc_dummy if filter_1 == 1 , percent(vacc_dummy) stack asyvars recast(bar) ///
		fxsize(25)									///
		title("Free grilled" "sausage")	///
		ytitle("")									///
		l1title("", )								///
		ylabel(,)					///
		blabel(bar, pos(center) format(%2.1f) )	///
		var1opts(gap(*0) label(labsize(large))) ///
			bar(1, color(538t) fintensity(inten90)) ///
			bar(2, color(538t) fintensity(inten50)) ///
			bar(3, color(gray) fintensity(inten40)) ///
			bar(4, color(538o) fintensity(inten50)) ///
			bar(5, color(538o) fintensity(inten90)) ///
		legend(rows(1) stack size(small) ///
		symplacement(center) pos(6) ///
		label(1 "not at all meaningful") 	///
		label(2 "rather not meaningful")	///
		label(3 "neutral") 	///
		label(4 "rather meaningful") 	///
		label(5 "very meaningful")) 	
		graph save "dump/fig_02b.gph", replace
		
	catplot incentive_lottery vacc_dummy if filter_1 == 1 , percent(vacc_dummy) stack asyvars recast(bar) ///
		fxsize(25)									///
		title("Vaccination lottery" "(chance to win big)")	///
		ytitle("")									///
		l1title("", )								///
		ylabel(,)					///
		blabel(bar, pos(center) format(%2.1f) )	///
		var1opts(gap(*0) label(labsize(large))) ///
			bar(1, color(538t) fintensity(inten90)) ///
			bar(2, color(538t) fintensity(inten50)) ///
			bar(3, color(gray) fintensity(inten40)) ///
			bar(4, color(538o) fintensity(inten50)) ///
			bar(5, color(538o) fintensity(inten90)) ///
		legend(rows(1) stack size(small) ///
		symplacement(center) pos(6) ///
		label(1 "not at all meaningful") 	///
		label(2 "rather not meaningful")	///
		label(3 "neutral") 	///
		label(4 "rather meaningful") 	///
		label(5 "very meaningful")) 	
		graph save "dump/fig_02c.gph", replace
		
	catplot incentive_money vacc_dummy if filter_1 == 1 , percent(vacc_dummy) stack asyvars recast(bar) ///
		fxsize(25)									///
		title("Direct monetary" "compensation")	///
		ytitle("")									///
		l1title("", )								///
		ylabel(,)					///
		blabel(bar, pos(center) format(%2.1f) )	///
		var1opts(gap(*0) label(labsize(large))) ///
			bar(1, color(538t) fintensity(inten90)) ///
			bar(2, color(538t) fintensity(inten50)) ///
			bar(3, color(gray) fintensity(inten40)) ///
			bar(4, color(538o) fintensity(inten50)) ///
			bar(5, color(538o) fintensity(inten90)) ///
		legend(rows(1) stack size(small) ///
		symplacement(center) pos(6) ///
		label(1 "not at all meaningful") 	///
		label(2 "rather not meaningful")	///
		label(3 "neutral") 	///
		label(4 "rather meaningful") 	///
		label(5 "very meaningful")) 	
		graph save "dump/fig_02d.gph", replace

		grc1leg 	"dump/fig_02a.gph" ///
					"dump/fig_02b.gph" ///
					"dump/fig_02c.gph" ///
					"dump/fig_02d.gph" , col(4) imargin (1 0 0 0) commonscheme title("Outcome I: How would you rate the following incentives" "to encourage people to get a COVID-19 vaccination?")

		graph save "dump/fig_02A.gph", replace			
		
		catplot   vacc_dummy money_incentive if filter_1 == 1 , percent(vacc_dummy) 				///
					asyvars recast(hbar) 												///
					legend(col(1) ring(0) position(5) bmargin(medium) size(vsmall) title("#SARS-CoV2 vacc." "doses received", size(vsmall)))		///
					fysize(40) 															///
					ylabel(,labsize(vsmall))											///
					ytitle("percent" ,size(small))										///
					var2opts( label(labsize(vsmall)))								///
					title("Outcome II: If there were a financial incentive for COVID-19 vaccination," "how much money would you consider appropriate?") 
		graph save "dump/fig_02B.gph", replace
		
		graph combine 	"dump/fig_02A.gph" ///
						"dump/fig_02B.gph", col(1) imargin (0 0 0 0) xsize(6) ysize(7)
						
						
		gr_edit 		.plotregion1.graph1.plotregion1.graph1.plotregion1.barlabels[2].DragBy -1.142280608047837 22.82679205889677

		gr_edit 		.plotregion1.graph1.plotregion1.graph1.plotregion1.barlabels[5].DragBy 0 22.82679205889677

		gr_edit 		.plotregion1.graph1.plotregion1.graph1.plotregion1.barlabels[4].DragBy -.380760202682603 -22.22608700471527

		gr_edit 		.plotregion1.graph1.plotregion1.graph1.plotregion1.barlabels[2].DragBy .7615204053653174 0

		gr_edit 		.plotregion1.graph1.plotregion1.graph2.plotregion1.barlabels[5].DragBy .1903801034086672 22.82679205889689

		gr_edit 		.plotregion1.graph1.plotregion1.graph2.plotregion1.barlabels[4].DragBy 0 -21.62538195053376

		gr_edit 		.plotregion1.graph1.plotregion1.graph2.plotregion1.barlabels[2].DragBy -.3807602068173899 22.22608700471528

		gr_edit 		.plotregion1.graph1.plotregion1.graph3.plotregion1.barlabels[5].DragBy 0 22.82679205889698

		gr_edit 		.plotregion1.graph1.plotregion1.graph3.plotregion1.barlabels[4].DragBy -.3807602075904133 -21.6253819505339

		gr_edit 		.plotregion1.graph1.plotregion1.graph3.plotregion1.barlabels[2].DragBy -.1903801037951928 22.82679205889689

		gr_edit 		.plotregion1.graph1.plotregion1.graph4.plotregion1.barlabels[5].DragBy -.1903801013413433 22.82679205889689

		gr_edit 		.plotregion1.graph1.plotregion1.graph4.plotregion1.barlabels[4].DragBy -.1903801013412876 -22.22608700471528

		gr_edit 		.plotregion1.graph1.plotregion1.graph4.plotregion1.barlabels[2].DragBy 0 22.22608700471528
		
		gr_edit			.AddTextBox added_text editor 50.83140785129991 -1.928044062301026
		gr_edit			.added_text_new = 1
		gr_edit			.added_text_rec = 1
		gr_edit			.added_text[1].style.editstyle  angle(default) size(medsmall) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
		gr_edit			.added_text[1].text = {}
		gr_edit			.added_text[1].text.Arrpush `"# SARS-CoV-2"'
		gr_edit			.added_text[1].text.Arrpush `"vaccination doses"'
		gr_edit			.added_text[1].style.editstyle size(tiny) editcopy
		gr_edit			.added_text[1].DragBy -1.087823750795866 1.813039584659828
		gr_edit			.added_text[1].DragBy -.1813039584659832 -.362607916931964
		gr_edit			.AddTextBox added_text editor 111.0709682260775 -.6897204236523002
		gr_edit			.added_text_new = 2
		gr_edit			.added_text_rec = 2
		gr_edit			.added_text[2].style.editstyle  angle(default) size(medlarge) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
		gr_edit			.added_text[2].text = {}
		gr_edit			.added_text[2].text.Arrpush a
		gr_edit			.AddTextBox added_text editor 38.3258352986633 -.8155770550146082
		gr_edit			.added_text_new = 3
		gr_edit			.added_text_rec = 3
		gr_edit			.added_text[3].style.editstyle  angle(default) size(medlarge) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
		gr_edit			.added_text[3].text = {}
		gr_edit			.added_text[3].text.Arrpush b
		gr_edit			.added_text[3].DragBy -.7551397881738351 .7551397881738504

		
/*
		gr_edit			.plotregion1.graph2.plotregion1.AddTextBox added_text editor 49.06024379708616 87.18130972469858
		gr_edit			.plotregion1.graph2.plotregion1.added_text_new = 1
		gr_edit			.plotregion1.graph2.plotregion1.added_text_rec = 1
		gr_edit			.plotregion1.graph2.plotregion1.added_text[1].style.editstyle  angle(default) size(medsmall) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
		gr_edit			.plotregion1.graph2.plotregion1.added_text[1].style.editstyle size(tiny) editcopy
		gr_edit			.plotregion1.graph2.plotregion1.added_text[1].text = {}
		gr_edit			.plotregion1.graph2.plotregion1.added_text[1].text.Arrpush `"# SARS-CoV-2"'
		gr_edit			.plotregion1.graph2.plotregion1.added_text[1].text.Arrpush `"vaccination doses"'
		gr_edit			.plotregion1.graph2.plotregion1.added_text[1].style.editstyle size(vsmall) editcopy
		gr_edit			.plotregion1.graph2.plotregion1.added_text[1].DragBy -11.30276230457268 0
		gr_edit			.plotregion1.graph2.plotregion1.added_text[1].DragBy 3.477773016791587 0
*/

		graph export "dump\fig_02_allgroups.emf", replace
		graph export "dump\fig_02_allgroups.pdf", replace		
		
	use "dump\PPD_July_2022_labeled.dta", replace
	
	catplot incentive_football vacc_dummy if filter_1 == 1 & framing == 0 , percent(vacc_dummy) stack asyvars recast(bar) ///
		fxsize(25)									///
		title("Free ticket for" "football match")	///
		ytitle("")									///
		l1title("", )								///
		ylabel(,)					///
		blabel(bar, pos(center) format(%2.1f) )	///
		var1opts(gap(*0) label(labsize(large))) ///
			bar(1, color(538t) fintensity(inten90)) ///
			bar(2, color(538t) fintensity(inten50)) ///
			bar(3, color(gray) fintensity(inten40)) ///
			bar(4, color(538o) fintensity(inten50)) ///
			bar(5, color(538o) fintensity(inten90)) ///
		legend(rows(1) stack size(vsmall) 	///
		symplacement(center) pos(6) 		///
		label(1 "not at all meaningful") 	///
		label(2 "rather not meaningful")	///
		label(3 "neutral") 	///
		label(4 "rather meaningful") 	///
		label(5 "very meaningful")) 	///
		text(50 -42 "percent", orientation(vertical) size(medsmall) )
		graph save "dump/fig_02a.gph", replace

	catplot incentive_sausage vacc_dummy if filter_1 == 1 & framing == 0 , percent(vacc_dummy) stack asyvars recast(bar) ///
		fxsize(25)									///
		title("Free grilled" "sausage")	///
		ytitle("")									///
		l1title("", )								///
		ylabel(,)					///
		blabel(bar, pos(center) format(%2.1f) )	///
		var1opts(gap(*0) label(labsize(large))) ///
			bar(1, color(538t) fintensity(inten90)) ///
			bar(2, color(538t) fintensity(inten50)) ///
			bar(3, color(gray) fintensity(inten40)) ///
			bar(4, color(538o) fintensity(inten50)) ///
			bar(5, color(538o) fintensity(inten90)) ///
		legend(rows(1) stack size(small) ///
		symplacement(center) pos(6) ///
		label(1 "not at all meaningful") 	///
		label(2 "rather not meaningful")	///
		label(3 "neutral") 	///
		label(4 "rather meaningful") 	///
		label(5 "very meaningful")) 	
		graph save "dump/fig_02b.gph", replace
		
	catplot incentive_lottery vacc_dummy if filter_1 == 1 & framing == 0 , percent(vacc_dummy) stack asyvars recast(bar) ///
		fxsize(25)									///
		title("Vaccination lottery" "(chance to win big)")	///
		ytitle("")									///
		l1title("", )								///
		ylabel(,)					///
		blabel(bar, pos(center) format(%2.1f) )	///
		var1opts(gap(*0) label(labsize(large))) ///
			bar(1, color(538t) fintensity(inten90)) ///
			bar(2, color(538t) fintensity(inten50)) ///
			bar(3, color(gray) fintensity(inten40)) ///
			bar(4, color(538o) fintensity(inten50)) ///
			bar(5, color(538o) fintensity(inten90)) ///
		legend(rows(1) stack size(small) ///
		symplacement(center) pos(6) ///
		label(1 "not at all meaningful") 	///
		label(2 "rather not meaningful")	///
		label(3 "neutral") 	///
		label(4 "rather meaningful") 	///
		label(5 "very meaningful")) 	
		graph save "dump/fig_02c.gph", replace
		
	catplot incentive_money vacc_dummy if filter_1 == 1 & framing == 0, percent(vacc_dummy) stack asyvars recast(bar) ///
		fxsize(25)									///
		title("Direct monetary" "compensation")	///
		ytitle("")									///
		l1title("", )								///
		ylabel(,)					///
		blabel(bar, pos(center) format(%2.1f) )	///
		var1opts(gap(*0) label(labsize(large))) ///
			bar(1, color(538t) fintensity(inten90)) ///
			bar(2, color(538t) fintensity(inten50)) ///
			bar(3, color(gray) fintensity(inten40)) ///
			bar(4, color(538o) fintensity(inten50)) ///
			bar(5, color(538o) fintensity(inten90)) ///
		legend(rows(1) stack size(small) ///
		symplacement(center) pos(6) ///
		label(1 "not at all meaningful") 	///
		label(2 "rather not meaningful")	///
		label(3 "neutral") 	///
		label(4 "rather meaningful") 	///
		label(5 "very meaningful")) 	
		graph save "dump/fig_02d.gph", replace

		grc1leg 	"dump/fig_02a.gph" ///
					"dump/fig_02b.gph" ///
					"dump/fig_02c.gph" ///
					"dump/fig_02d.gph" , col(4) imargin (1 0 0 0) commonscheme title("Outcome I: How would you rate the following incentives" "to encourage people to get a COVID-19 vaccination?")

		graph save "dump/fig_02A.gph", replace			
		
		catplot   vacc_dummy money_incentive if filter_1 == 1 & framing == 0, percent(vacc_dummy) 				///
					asyvars recast(hbar) 												///
					legend(col(1) ring(0) position(5) bmargin(medium) size(vsmall) title("#SARS-CoV2 vacc." "doses received", size(vsmall)))		///
					fysize(40) 															///
					ylabel(,labsize(vsmall))											///
					ytitle("percent" ,size(small))										///
					var2opts( label(labsize(vsmall)))								///
					title("Outcome II: If there were a financial incentive for COVID-19 vaccination," "how much money would you consider appropriate?") 
		graph save "dump/fig_02B.gph", replace
		
		graph combine 	"dump/fig_02A.gph" ///
						"dump/fig_02B.gph", col(1) imargin (0 0 0 0) xsize(6) ysize(7)
						
						
		gr_edit 		.plotregion1.graph1.plotregion1.graph1.plotregion1.barlabels[2].DragBy -1.142280608047837 22.82679205889677

		gr_edit 		.plotregion1.graph1.plotregion1.graph1.plotregion1.barlabels[5].DragBy 0 22.82679205889677

		gr_edit 		.plotregion1.graph1.plotregion1.graph1.plotregion1.barlabels[4].DragBy -.380760202682603 -22.22608700471527

		gr_edit 		.plotregion1.graph1.plotregion1.graph1.plotregion1.barlabels[2].DragBy .7615204053653174 0

		gr_edit 		.plotregion1.graph1.plotregion1.graph2.plotregion1.barlabels[5].DragBy .1903801034086672 22.82679205889689

		gr_edit 		.plotregion1.graph1.plotregion1.graph2.plotregion1.barlabels[4].DragBy 0 -21.62538195053376

		gr_edit 		.plotregion1.graph1.plotregion1.graph2.plotregion1.barlabels[2].DragBy -.3807602068173899 22.22608700471528

		gr_edit 		.plotregion1.graph1.plotregion1.graph3.plotregion1.barlabels[5].DragBy 0 22.82679205889698

		gr_edit 		.plotregion1.graph1.plotregion1.graph3.plotregion1.barlabels[4].DragBy -.3807602075904133 -21.6253819505339

		gr_edit 		.plotregion1.graph1.plotregion1.graph3.plotregion1.barlabels[2].DragBy -.1903801037951928 22.82679205889689

		gr_edit 		.plotregion1.graph1.plotregion1.graph4.plotregion1.barlabels[5].DragBy -.1903801013413433 22.82679205889689

		gr_edit 		.plotregion1.graph1.plotregion1.graph4.plotregion1.barlabels[4].DragBy -.1903801013412876 -22.22608700471528

		gr_edit 		.plotregion1.graph1.plotregion1.graph4.plotregion1.barlabels[2].DragBy 0 22.22608700471528
		
		gr_edit			.AddTextBox added_text editor 50.83140785129991 -1.928044062301026
		gr_edit			.added_text_new = 1
		gr_edit			.added_text_rec = 1
		gr_edit			.added_text[1].style.editstyle  angle(default) size(medsmall) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
		gr_edit			.added_text[1].text = {}
		gr_edit			.added_text[1].text.Arrpush `"# SARS-CoV-2"'
		gr_edit			.added_text[1].text.Arrpush `"vaccination doses"'
		gr_edit			.added_text[1].style.editstyle size(tiny) editcopy
		gr_edit			.added_text[1].DragBy -1.087823750795866 1.813039584659828
		gr_edit			.added_text[1].DragBy -.1813039584659832 -.362607916931964
		gr_edit			.AddTextBox added_text editor 111.0709682260775 -.6897204236523002
		gr_edit			.added_text_new = 2
		gr_edit			.added_text_rec = 2
		gr_edit			.added_text[2].style.editstyle  angle(default) size(medlarge) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
		gr_edit			.added_text[2].text = {}
		gr_edit			.added_text[2].text.Arrpush a
		gr_edit			.AddTextBox added_text editor 38.3258352986633 -.8155770550146082
		gr_edit			.added_text_new = 3
		gr_edit			.added_text_rec = 3
		gr_edit			.added_text[3].style.editstyle  angle(default) size(medlarge) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
		gr_edit			.added_text[3].text = {}
		gr_edit			.added_text[3].text.Arrpush b
		gr_edit			.added_text[3].DragBy -.7551397881738351 .7551397881738504

		
/*
		gr_edit			.plotregion1.graph2.plotregion1.AddTextBox added_text editor 49.06024379708616 87.18130972469858
		gr_edit			.plotregion1.graph2.plotregion1.added_text_new = 1
		gr_edit			.plotregion1.graph2.plotregion1.added_text_rec = 1
		gr_edit			.plotregion1.graph2.plotregion1.added_text[1].style.editstyle  angle(default) size(medsmall) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
		gr_edit			.plotregion1.graph2.plotregion1.added_text[1].style.editstyle size(tiny) editcopy
		gr_edit			.plotregion1.graph2.plotregion1.added_text[1].text = {}
		gr_edit			.plotregion1.graph2.plotregion1.added_text[1].text.Arrpush `"# SARS-CoV-2"'
		gr_edit			.plotregion1.graph2.plotregion1.added_text[1].text.Arrpush `"vaccination doses"'
		gr_edit			.plotregion1.graph2.plotregion1.added_text[1].style.editstyle size(vsmall) editcopy
		gr_edit			.plotregion1.graph2.plotregion1.added_text[1].DragBy -11.30276230457268 0
		gr_edit			.plotregion1.graph2.plotregion1.added_text[1].DragBy 3.477773016791587 0
*/

		graph export "fig\fig_02.emf", replace
		graph export "fig\fig_02.pdf", replace		
		
		
	************************************************************************	
	* Figure 03: Average money deemed acceptable as compensation for vaccination by the attitude towards money as an incentive for vaccination in general 
	************************************************************************
	use "dump\PPD_July_2022_labeled.dta", replace

		graph hbar (mean) money_incentive_metric if filter_1 ==1, over(incentive_money) ytitle("EUR", size(huge)) xsize(6) ysize(1) scale(1.5) 
			gr_edit scaleaxis.reset_rule 20, tickset(major) ruletype(suggest) 
			gr_edit grpaxis.major.num_rule_ticks = 0
			gr_edit grpaxis.edit_tick 4 70.6683 `"rather not meaningful"', tickset(major)
			gr_edit grpaxis.major.num_rule_ticks = 0
			gr_edit grpaxis.edit_tick 3 50 `"neutral"', tickset(major)
			gr_edit grpaxis.major.num_rule_ticks = 0
			gr_edit grpaxis.edit_tick 2 29.3317 `"rather meaningful"', tickset(major)
			gr_edit grpaxis.style.editstyle majorstyle(tickstyle(textstyle(size(huge)))) editcopy
			gr_edit scaleaxis.style.editstyle majorstyle(tickstyle(textstyle(size(huge)))) editcopy
		graph export "dump\fig_03_allgroups.emf", replace	
		graph export "dump\fig_03_allgroups.pdf", replace		
		
		use "dump\PPD_July_2022_labeled.dta", replace

		graph hbar (mean) money_incentive_metric if filter_1 ==1 & framing == 0, over(incentive_money) ytitle("EUR", size(huge)) xsize(6) ysize(1) scale(1.5) 
			gr_edit scaleaxis.reset_rule 20, tickset(major) ruletype(suggest) 
			gr_edit grpaxis.major.num_rule_ticks = 0
			gr_edit grpaxis.edit_tick 4 70.6683 `"rather not meaningful"', tickset(major)
			gr_edit grpaxis.major.num_rule_ticks = 0
			gr_edit grpaxis.edit_tick 3 50 `"neutral"', tickset(major)
			gr_edit grpaxis.major.num_rule_ticks = 0
			gr_edit grpaxis.edit_tick 2 29.3317 `"rather meaningful"', tickset(major)
			gr_edit grpaxis.style.editstyle majorstyle(tickstyle(textstyle(size(huge)))) editcopy
			gr_edit scaleaxis.style.editstyle majorstyle(tickstyle(textstyle(size(huge)))) editcopy
		graph export "fig\fig_03.emf", replace	
		graph export "fig\fig_03.pdf", replace		
		
	
		
	********************************************************************************************************************************************************************	
	* Fig 04: Differences in the means of the dependent variables between the framing groups + results of Dunnett’s test
	********************************************************************************************************************************************************************
		* Dunnet's test
		* Read data 
	use "dump\PPD_July_2022_labeled.dta", replace
	drop if filter_1 == 0
	bysort framing: sum(incentive_football)
	pwmean incentive_football , over(framing)  mcompare(dunnett) effects
	** --> no significant differences between experimental groups and control group
	
	bysort framing: sum(incentive_sausage)
	pwmean incentive_sausage , over(framing)  mcompare(dunnett) effects
	** --> no significant differences between experimental groups and control group
	
	bysort framing: sum(incentive_lottery)	
	pwmean incentive_lottery , over(framing)  mcompare(dunnett) effects
	** --> Long Covid Frame 2: p < 0.05; Economic Frame 2: p < 0.01
	
	bysort framing: sum(incentive_money)	
	pwmean incentive_money , over(framing)  mcompare(dunnett) effects
	** --> Long Covid Frame 2: p < 0.05; Economic Frame 1: p < 0.1
	
	
		* Read data 
	use "dump\PPD_July_2022_labeled.dta", replace
	
	* save means in new dataset
	collapse incentive_* if filter_1 == 1, by(framing)
	* transpose dataset
	xpose , clear varname
	drop if _varname == "framing"
	label variable v1 "Control"
	label variable v2 "Long COVID frame 1"
	label variable v3 "Long COVID-frame 2"
	label variable v4 "Economic frame 1"
	label variable v5 "Economic frame 2"

	gen incentive = .
	forvalues i = 1/4 {
	replace incentive = `i' in `i'
	}
	label define incentive 1 "Football ticket" 2 "Sausage" 3 "Lottery" 4 "Money", modify
	label value incentive incentive

	graph dot (asis) v*, over(incentive)  							///
					ylabel(.15(.05).35) exclude0 					/// limit range of x-axis
					legend(pos(6) rows(1)) 							/// legend in one row below
					fysize(60)										///
					ytitle("Average meaningfulness (scale 0 to 1)") 		///
					marker(1, msize(3)) 							///
					marker(2, msize(3)) 							///
					marker(3, msize(3)) 							///
					marker(4, msize(3)) 							///
					marker(5, msize(3)) scale(1.2) 					///
					text (0.201 90  "*") 						///
					text (0.225 15.3 "**" )						///
					text (0.22 15.3 "*") 							///
					text (0.319 40.4 "**")							///
					text (0.326 40.4 "***")						///
					title("Outcome I: Attitudes towards incentives") 
			*		caption("Stars indicate significant differences to control group calculated via Dunnett's test. * p < 0.1; ** p < 0.05; *** p < 0.01.", size(small)) ///

	graph save "dump\fig_A04a.gph", replace	
	
	* Read data 
	use "dump\PPD_July_2022_labeled.dta", replace
	drop if filter_1 == 0
	bysort framing: sum(money_incentive_metric)	
	pwmean money_incentive_metric , over(framing)  mcompare(dunnett) effects
	** --> Long Covid Frame 2: p < 0.1

	* save means in new dataset
	collapse money_incentive_metric if filter_1 == 1, by(framing)
	* transpose dataset
	xpose , clear varname
	drop if _varname == "framing"
	label variable v1 "Control"
	label variable v2 "Long COVID frame 1"
	label variable v3 "Long COVID-frame 2"
	label variable v4 "Economic frame 1"
	label variable v5 "Economic frame 2"

	gen money = .
	replace money = 1 in 1

	label define money 1 "How much money?" , modify
	label value money money


	graph dot (asis) v*, over(money)  							///
					ylabel(22(2)32) exclude0 					/// limit range of x-axis
					legend(pos(6) rows(1)) 						/// legend in one row below
					fysize(30) 							///
					ytitle("Compensation for getting COVID-19 vaccine (group means in EUR)") 	///
					marker(1, msize(3)) 							///
					marker(2, msize(3)) 							///
					marker(3, msize(3)) 							///
					marker(4, msize(3)) 							///
					marker(5, msize(3)) scale(1.2) 					///
					text (30.57 66  "*") 						///
					///caption("Stars indicate significant differences to control group calculated via Dunnett's test. * p < 0.1; ** p < 0.05; *** p < 0.01.", size(mediumsmall)) ///
					title("Outcome II: How much money should be paid as incentive to get vaccinated?")

	graph save "dump\fig_A04b.gph", replace	

	grc1leg "dump\fig_A04a.gph" "dump\fig_A04b.gph", row(2) scale(.75) imargin(0 0 20 0) 
	gr_edit .plotregion1.graph2.caption.DragBy -17.73589921627527 -8.513064051108556
	gr_edit .legend.DragBy 9.739462523929522 -.1064959168722647
	gr_edit .plotregion1.graph2.caption.DragBy 1.023224954785111 0
	gr_edit .plotregion1.graph1.grpaxis.style.editstyle majorstyle(tickstyle(textstyle(size(medium)))) editcopy
	gr_edit .plotregion1.graph2.grpaxis.style.editstyle majorstyle(tickstyle(textstyle(size(medium)))) editcopy
	gr_edit .plotregion1.graph2.caption.DragBy 5.974283761888121 -6.521841034200942
	gr_edit .plotregion1.graph1.title.style.editstyle size(medium) editcopy
	gr_edit .plotregion1.graph2.title.style.editstyle size(medium) editcopy

	graph export "fig\fig_04.emf", replace			
	graph export "fig\fig_04.pdf", replace
	
	
	
	
		

	********************************************************************************************************************************************************************	
	* Fig 05: Analysis I - Meaningfulness of vaccination incentives predicted based on multiple OLS-model
	********************************************************************************************************************************************************************		
	
	use "dump\PPD_July_2022_labeled.dta", replace
	
/*
	reg incentive_football  													/// Analysis I DV 
			i.framing															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	margins, at(framing=(0(1)4)  ) saving(football, replace)
	
	reg incentive_sausage 													/// Analysis I DV 
			i.framing															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	margins, at(framing=(0(1)4)  ) saving(sausage, replace)

	
	reg incentive_lottery  													/// Analysis I DV 
			i.framing															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	margins, at(framing=(0(1)4)  ) saving(lottery, replace)

	
	reg incentive_money 													/// Analysis I DV 
			i.framing															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 
	margins, at(framing=(0(1)4)  ) saving(money, replace)
*/

	graph combine "dump/fig_football_ols.gph" "dump/fig_sausage_ols.gph" "dump/fig_lottery_ols.gph" "dump/fig_money_ols.gph", row(1) xsize(8) ysize(3) scale(1.5)
	gr_edit .plotregion1.graph1.xaxis1.major.num_rule_ticks = 0
	gr_edit .plotregion1.graph1.xaxis1.edit_tick 1 0 `"Control group"', custom tickset(major) editstyle(tickstyle(textstyle(color(reddish))) )
	gr_edit .plotregion1.graph2.xaxis1.major.num_rule_ticks = 0
	gr_edit .plotregion1.graph2.xaxis1.edit_tick 1 0 `"Control group"', custom tickset(major) editstyle(tickstyle(textstyle(color(reddish))) )
	gr_edit .plotregion1.graph3.xaxis1.major.num_rule_ticks = 0
	gr_edit .plotregion1.graph3.xaxis1.edit_tick 1 0 `"Control group"', custom tickset(major) editstyle(tickstyle(textstyle(color(reddish))) )
	gr_edit .plotregion1.graph4.xaxis1.major.num_rule_ticks = 0
	gr_edit .plotregion1.graph4.xaxis1.edit_tick 1 0 `"Control group"', custom tickset(major) editstyle(tickstyle(textstyle(color(reddish))) )
	gr_edit .plotregion1.graph1.plotregion1.AddMarker added_markers editor .0054740657037635 .2233700283787271
	gr_edit .plotregion1.graph1.plotregion1.added_markers_new = 1
	gr_edit .plotregion1.graph1.plotregion1.added_markers_rec = 1
	gr_edit .plotregion1.graph1.plotregion1.added_markers[1].style.editstyle  marker( symbol(square) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(reddish) size(medsmall) angle(stdarrow) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(zero) backangle(stdarrow) backsymangle(zero)) line( width(thin) color(black) pattern(solid) align(inside)) area( linestyle( width(thin) color(white) pattern(solid) align(inside)) shadestyle( color(white) intensity(inten100) fill(pattern10))) label( textstyle( horizontal(center) vertical(middle) angle(default) size(small) color(black) position() margin(zero) linestyle( width(thin) color(black) pattern(solid) align(inside))) position(6) textgap(half_tiny) horizontal(default) vertical(default)) dots( symbol(circle) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(black) size(vtiny) angle(horizontal) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(small) backangle(horizontal) backsymangle(zero)) connect(direct) connect_missings(yes) editcopy
	gr_edit .plotregion1.graph2.plotregion1.AddMarker added_markers editor -.009039258484574 .2064617743808659
	gr_edit .plotregion1.graph2.plotregion1.added_markers_new = 1
	gr_edit .plotregion1.graph2.plotregion1.added_markers_rec = 1
	gr_edit .plotregion1.graph2.plotregion1.added_markers[1].style.editstyle  marker( symbol(square) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(reddish) size(medsmall) angle(stdarrow) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(zero) backangle(stdarrow) backsymangle(zero)) line( width(thin) color(black) pattern(solid) align(inside)) area( linestyle( width(thin) color(white) pattern(solid) align(inside)) shadestyle( color(white) intensity(inten100) fill(pattern10))) label( textstyle( horizontal(center) vertical(middle) angle(default) size(small) color(black) position() margin(zero) linestyle( width(thin) color(black) pattern(solid) align(inside))) position(6) textgap(half_tiny) horizontal(default) vertical(default)) dots( symbol(circle) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(black) size(vtiny) angle(horizontal) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(small) backangle(horizontal) backsymangle(zero)) connect(direct) connect_missings(yes) editcopy
	gr_edit .plotregion1.graph3.plotregion1.AddMarker added_markers editor .0054740657037635 .2773912235717336
	gr_edit .plotregion1.graph3.plotregion1.added_markers_new = 1
	gr_edit .plotregion1.graph3.plotregion1.added_markers_rec = 1
	gr_edit .plotregion1.graph3.plotregion1.added_markers[1].style.editstyle  marker( symbol(square) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(reddish) size(medsmall) angle(stdarrow) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(zero) backangle(stdarrow) backsymangle(zero)) line( width(thin) color(black) pattern(solid) align(inside)) area( linestyle( width(thin) color(white) pattern(solid) align(inside)) shadestyle( color(white) intensity(inten100) fill(pattern10))) label( textstyle( horizontal(center) vertical(middle) angle(default) size(small) color(black) position() margin(zero) linestyle( width(thin) color(black) pattern(solid) align(inside))) position(6) textgap(half_tiny) horizontal(default) vertical(default)) dots( symbol(circle) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(black) size(vtiny) angle(horizontal) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(small) backangle(horizontal) backsymangle(zero)) connect(direct) connect_missings(yes) editcopy
	gr_edit .plotregion1.graph4.plotregion1.AddMarker added_markers editor -.009039258484574 .1905955197160464
	gr_edit .plotregion1.graph4.plotregion1.added_markers_new = 1
	gr_edit .plotregion1.graph4.plotregion1.added_markers_rec = 1
	gr_edit .plotregion1.graph4.plotregion1.added_markers[1].style.editstyle  marker( symbol(square) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(reddish) size(medsmall) angle(stdarrow) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(zero) backangle(stdarrow) backsymangle(zero)) line( width(thin) color(black) pattern(solid) align(inside)) area( linestyle( width(thin) color(white) pattern(solid) align(inside)) shadestyle( color(white) intensity(inten100) fill(pattern10))) label( textstyle( horizontal(center) vertical(middle) angle(default) size(small) color(black) position() margin(zero) linestyle( width(thin) color(black) pattern(solid) align(inside))) position(6) textgap(half_tiny) horizontal(default) vertical(default)) dots( symbol(circle) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(black) size(vtiny) angle(horizontal) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(small) backangle(horizontal) backsymangle(zero)) connect(direct) connect_missings(yes) editcopy

	graph export "fig\fig_05.emf", replace			
	graph export "fig\fig_05.pdf", replace	

	*** OLD Version
// 	combomarginsplot12 football sausage lottery money, labels("Football" "Sausage" "Lottery" "Money")  recast(scatter) savefile(combined_margins, replace)
//	
// 	combomarginsplot12 football, recast(scatter) ytitle("Meaningfulness of vaccine incentives" "(0 = not at all; 1 = very meaningful)")  title("Football")
//	
// 	combomarginsplot12 football sausage lottery money, labels("Football" "Sausage" "Lottery" "Money") horiz recast(scatter) savefile(combined_margins, replace)
// 	mplotoffset using combined_margins,  recast(scatter) xtitle("") ytitle("Meaningfulness of vaccine incentives" "(0 = not at all; 1 = very meaningful)")  title("") legend(pos(6) row(1))
// 		gr_edit .plotregion1.AddLine added_lines editor .0234340238741471 .2769536618514004 4.126906308084747 .2769536618514004
// 		gr_edit .plotregion1.added_lines_new = 1
// 		gr_edit .plotregion1.added_lines_rec = 1
// 		gr_edit .plotregion1.added_lines[1].style.editstyle  linestyle( width(thin) color(green) pattern(longdash) align(inside)) headstyle( symbol(circle) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(black) size(medium) angle(stdarrow) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(zero) backangle(stdarrow) backsymangle(zero)) headpos(neither) editcopy
//
//
// 		gr_edit .plotregion1.AddLine added_lines editor -.0733459828289329 .2238233058888829 4.126906308084743 .2238233058888829
// 		gr_edit .plotregion1.added_lines_new = 2
// 		gr_edit .plotregion1.added_lines_rec = 2
// 		gr_edit .plotregion1.added_lines[2].style.editstyle  linestyle( width(thin) color(green) pattern(longdash) align(inside)) headstyle( symbol(circle) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(black) size(medium) angle(stdarrow) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(zero) backangle(stdarrow) backsymangle(zero)) headpos(neither) editcopy
//
//
// 		gr_edit .plotregion1.AddLine added_lines editor -.0230203793433313 .2062143879127341 4.111421507012252 .2062143879127341
// 		gr_edit .plotregion1.added_lines_new = 3
// 		gr_edit .plotregion1.added_lines_rec = 3
// 		gr_edit .plotregion1.added_lines[3].style.editstyle  linestyle( width(thin) color(green) pattern(longdash) align(inside)) headstyle( symbol(circle) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(black) size(medium) angle(stdarrow) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(zero) backangle(stdarrow) backsymangle(zero)) headpos(neither) editcopy
//
//
// 		gr_edit .plotregion1.AddLine added_lines editor .0737596273597488 .1904270821410146 4.107550306744128 .1904270821410146
// 		gr_edit .plotregion1.added_lines_new = 4
// 		gr_edit .plotregion1.added_lines_rec = 4
// 		gr_edit .plotregion1.added_lines[4].style.editstyle  linestyle( width(thin) color(green) pattern(longdash) align(inside)) headstyle( symbol(circle) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(black) size(medium) angle(stdarrow) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(zero) backangle(stdarrow) backsymangle(zero)) headpos(neither) editcopy
// 		// edits
//
// 		gr_edit .plotregion1.added_lines[1].style.editstyle linestyle(color(sky)) editcopy
// 		gr_edit .plotregion1.added_lines[1].style.editstyle headstyle(fillcolor(sky)) editcopy
// 		gr_edit .plotregion1.added_lines[1].style.editstyle headstyle(linestyle(color(sky))) editcopy
// 		// line[1] edits
//
// 		gr_edit .plotregion1.added_lines[2].style.editstyle linestyle(color(black)) editcopy
// 		gr_edit .plotregion1.added_lines[2].style.editstyle headstyle(fillcolor(black)) editcopy
// 		gr_edit .plotregion1.added_lines[2].style.editstyle headstyle(linestyle(color(black))) editcopy
// 		// line[2] edits
//
// 		gr_edit .plotregion1.added_lines[3].style.editstyle linestyle(color(gs10)) editcopy
// 		gr_edit .plotregion1.added_lines[3].style.editstyle headstyle(fillcolor(gs10)) editcopy
// 		gr_edit .plotregion1.added_lines[3].style.editstyle headstyle(linestyle(color(gs10))) editcopy
// 		// line[3] edits
//
// 		gr_edit .plotregion1.added_lines[1].drag_point = (.)
// 		gr_edit .plotregion1.added_lines[1].DragBy .0003036020340715 .0038712002681231
// 		// line[1] edits
//
// 		gr_edit .plotregion1.added_lines[2].drag_point = (.)
// 		gr_edit .plotregion1.added_lines[2].DragBy -.0006072040681431 0	
// 	graph export "fig\fig_05.emf", replace			
// 	graph export "fig\fig_05.pdf", replace	
	
	********************************************************************************************************************************************************************	
	* Fig 06: OLS regression explaining attitudes towards money and lottery as vaccination incentives (b-coefficients + 95% CI)
	********************************************************************************************************************************************************************		
	
	coefplot m1_ols_lottery, bylabel(Lottery)   || m1_ols_money, bylabel(Money)   || , xline(0) drop(_cons *.religion *.bundesland *.hh_income   )  	///
							headings(1.framing = "{it: Frames} (reference: control group)" 							///
							1.vacc_04 = "{it: Number of COVID-19 vaccinations received} (reference none)"			///
							school_system_3 = "{it: Attitudes towards medical and alternative belief systems}" 		///
							left_right = "{it: Political ideology}" 												///
							2.gender = "{it: Gender} (reference: male)" 											///
							2.age_group_raking = "{it: Age group} (reference: 18-30)" 								///
							2.sonntagsfrage = "{it: Voting intention} (reference: CDU/CSU)" 						///
							2.edu_3 = "{it: Education-Level} (reference: low)"										///
							1.bundesland = "{it: Bundesland} (reference: NRW)" 										///
							2.religion = "{it: Religion} (reference: nondenominational)"   							///
							1.vacc_04 = "{it: Vaccination status} (reference: not vaccinated)"						///
							vacc_decision_1 = "{it: The decision process to get or not get vaccinated ...}"			///
							threat_covid = "{it: Attitudes towards COVID-19}"										///
							2.hh_income = "{it: Household income} (reference: I find it very difficult to make ends meet)"																			///
							solidarity_index = "{it: Psychological personality}"									///
							spaltungsindex = " "																	///
							health_01 = " ")																		///
							xsize(10) ysize(8) 																		///
							xtitle("b-coefficient + 95% CI from multiple OLS regression")
							gr_edit .gmetric_mult = .6
							gr_edit .b1title.DragBy .3370319001386961 14.15214979195563
							gr_edit .b1title.style.editstyle size(vsmall) editcopy
							gr_edit .b1title.style.editstyle size(small) editcopy
							gr_edit .b1title.DragBy .5055478502080437 1.440638002773924

	graph export "fig\fig_06.emf", replace
	graph export "fig\fig_06.pdf", replace

	
	********************************************************************************************************************************************************************	
	* Fig 07: Money deemed acceptable as incentive for COVID-19 vaccination
	********************************************************************************************************************************************************************		
	
	use "dump\PPD_July_2022_labeled.dta", replace
	reg money_incentive_metric 													/// Analysis II DV 
			i.framing															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 

	margins , at(framing=(0(1)4)  )
	marginsplot, recast(scatter) title("") ytitle("") xtitle("Money to be paid as incentive for COVID-19 vaccination (in EUR)") horiz xsize(5) ysize(1.2) scale(2) level(95)
	graph export "fig\fig_07.emf", replace
	graph export "fig\fig_07.pdf", replace
	

	********************************************************************************************************************************************************************	
	* Fig 08: Predicted meaningfulness of vaccination incentives based on interaction models: frame x age/gender/vacc doses 
	********************************************************************************************************************************************************************	
	
	graph combine	"dump/fig_int_age_a.gph" "dump/fig_int_age_b.gph" "dump/fig_int_age_c.gph" "dump/fig_int_age_d.gph" , ycom xsize(10) ysize(6) rows(1) 
	graph save "dump\fig_08a.gph", replace
	
	graph combine	"dump/fig_int_gender_a.gph" "dump/fig_int_gender_b.gph" "dump/fig_int_gender_c.gph" "dump/fig_int_gender_d.gph" , ycom xsize(10) ysize(6) rows(1) 
	graph save "dump\fig_08b.gph", replace
	
	grc1leg2	"dump/fig_int_vacc_04_a.gph" "dump/fig_int_vacc_04_b.gph" "dump/fig_int_vacc_04_c.gph" "dump/fig_int_vacc_04_d.gph" , ycom xsize(10) ysize(6) rows(1) 
	graph save "dump\fig_08c.gph", replace
	

	
	///graph combine "dump\fig_08a.gph" "dump\fig_08b.gph" "dump\fig_08c.gph" "dump\fig_08d.gph", rows(4) xsize(6) ysize(8) scale(.9) 
	graph combine "dump\fig_08a.gph" "dump\fig_08b.gph" "dump\fig_08c.gph" , rows(3) xsize(10) ysize(8) scale(.9) 

	
	gr_edit .AddTextBox added_text editor 49.02551317614423 -1.548099861303745
	gr_edit .added_text_new = 1
	gr_edit .added_text_rec = 1
	gr_edit .added_text[1].style.editstyle  angle(default) size(vsmall) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
	gr_edit .added_text[1].style.editstyle size(vsmall) editcopy
	gr_edit .added_text[1]._set_orientation vertical
	gr_edit .added_text[1].text = {}
	gr_edit .added_text[1].text.Arrpush Meaningfulness of vaccine incentives (0 = not at all; 1 = very meaningful)
	gr_edit .plotregion1.graph3.legend.DragBy -2.253814147018027 -.4334257975034711
	
	gr_edit .AddTextBox added_text editor 65 60
	gr_edit .added_text_new = 2
	gr_edit .added_text_rec = 2
	gr_edit .added_text[2].style.editstyle  angle(default) size(vsmall) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
	gr_edit .added_text[2].style.editstyle size(vsmall) editcopy
	gr_edit .added_text[2]._set_orientation horizontal
	gr_edit .added_text[2].text = {}
	gr_edit .added_text[2].text.Arrpush Age
	
	gr_edit .AddTextBox added_text editor 33.5 58
	gr_edit .added_text_new = 3
	gr_edit .added_text_rec = 3
	gr_edit .added_text[3].style.editstyle  angle(default) size(vsmall) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
	gr_edit .added_text[3].style.editstyle size(vsmall) editcopy
	gr_edit .added_text[3]._set_orientation horizontal
	gr_edit .added_text[3].text = {}
	gr_edit .added_text[3].text.Arrpush Gender

	gr_edit .AddTextBox added_text editor 6 51
	gr_edit .added_text_new = 4
	gr_edit .added_text_rec = 4
	gr_edit .added_text[4].style.editstyle  angle(default) size(vsmall) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
	gr_edit .added_text[4].style.editstyle size(vsmall) editcopy
	gr_edit .added_text[4]._set_orientation horizontal
	gr_edit .added_text[4].text = {}
	gr_edit .added_text[4].text.Arrpush # SARS-CoV-2 vaccine doses	
	
	graph export "fig/fig_08.emf", replace
	graph export "fig/fig_08.pdf", replace	
	
******************************************************************
**																**
**		5)	Tables Main Article									**
**																**
******************************************************************	


******************************************************************
**																**
**		6)	Additional calculations								**
**																**
******************************************************************

	use "dump\PPD_July_2022_labeled.dta", replace
	* Mean money_incentive_metric all cases
	bysort gender: sum money_incentive_metric 
	tab money_incentive
	
	* overall meaningfulness of incentives
	tab incentive_football if filter_1 ==1
	tab incentive_sausage if filter_1 ==1
	tab incentive_lottery if filter_1 ==1
	tab incentive_money if filter_1 ==1
	
	* meaningfulness of incentives by vacc_status
	bysort vacc_dummy: tab incentive_lottery if filter_1 ==1
	bysort vacc_dummy: tab incentive_money if filter_1 ==1	
	
	* money_incentive by vacc_dummy
	tab  money_incentive vacc_dumm if filter_1 == 1, col nofreq
	
	* Mean money_incentive_metric only those who would pay something
	bysort vacc_dummy: sum money_incentive_metric if money_incentive != 1 & filter_1 == 1

	* Mean money_incentive_metric by incentive_money
	bysort incentive_money: sum money_incentive_metric if filter_1 ==1
	
	* Mean money_incentive_metric by framing for all
	bysort framing: sum money_incentive_metric if filter_1 ==1

	* Mean money_incentive_metric by framing for those who did not answer "nothing"
	bysort framing: sum money_incentive_metric if money_incentive != 1 & filter_1 == 1
	
	* 
	bysort framing: sum incentive_football if filter_1 ==1
	bysort framing: sum incentive_sausage if filter_1 ==1
	bysort framing: sum incentive_lottery if filter_1 ==1
	bysort framing: sum incentive_money if filter_1 ==1
	

	
	
******************************************************************
**																**
**		7)	Figures Online Annex								**
**																**
******************************************************************


	*****************************************************************
	* Fig S01: Overall acceptance of incentives by sex and age		
	*****************************************************************

	use "dump\PPD_July_2022_labeled.dta", replace
		
		catplot incentive_football gender if gender < 3 & filter_1 == 1 & framing == 0, percent(gender) stack asyvars recast(bar) ///
			fxsize(25)									///
			title("Free ticket for" "football match")	///
			ytitle("")									///
			l1title("", )								///
			ylabel(,)					///
			blabel(bar, pos(center) format(%2.1f) )	///
			var1opts(gap(*0) label(labsize(large))) ///
				bar(1, color(538t) fintensity(inten90)) ///
				bar(2, color(538t) fintensity(inten50)) ///
				bar(3, color(gray) fintensity(inten40)) ///
				bar(4, color(538o) fintensity(inten50)) ///
				bar(5, color(538o) fintensity(inten90)) ///
			legend(rows(1) stack size(vsmall) 	///
			symplacement(center) pos(6) 		///
			label(1 "not at all meaningful") 	///
			label(2 "rather not meaningful")	///
			label(3 "neutral") 	///
			label(4 "rather meaningful") 	///
			label(5 "very meaningful")) 	///
			text(50 -42 "percent", orientation(vertical) size(medsmall) )
			graph save "dump/fig_S01a.gph", replace

		catplot incentive_sausage gender if gender < 3 & filter_1 == 1 & framing == 0, percent(gender) stack asyvars recast(bar) ///
			fxsize(25)									///
			title("Free grilled" "sausage")	///
			ytitle("")									///
			l1title("", )								///
			ylabel(,)					///
			blabel(bar, pos(center) format(%2.1f) )	///
			var1opts(gap(*0) label(labsize(large))) ///
				bar(1, color(538t) fintensity(inten90)) ///
				bar(2, color(538t) fintensity(inten50)) ///
				bar(3, color(gray) fintensity(inten40)) ///
				bar(4, color(538o) fintensity(inten50)) ///
				bar(5, color(538o) fintensity(inten90)) ///
			legend(rows(1) stack size(small) ///
			symplacement(center) pos(6) ///
			label(1 "not at all meaningful") 	///
			label(2 "rather not meaningful")	///
			label(3 "neutral") 	///
			label(4 "rather meaningful") 	///
			label(5 "very meaningful")) 	
			graph save "dump/fig_S01b.gph", replace
			
		catplot incentive_lottery gender if gender < 3 & filter_1 == 1 & framing == 0, percent(gender) stack asyvars recast(bar) ///
			fxsize(25)									///
			title("Vaccination lottery" "(chance to win big)")	///
			ytitle("")									///
			l1title("", )								///
			ylabel(,)					///
			blabel(bar, pos(center) format(%2.1f) )	///
			var1opts(gap(*0) label(labsize(large))) ///
				bar(1, color(538t) fintensity(inten90)) ///
				bar(2, color(538t) fintensity(inten50)) ///
				bar(3, color(gray) fintensity(inten40)) ///
				bar(4, color(538o) fintensity(inten50)) ///
				bar(5, color(538o) fintensity(inten90)) ///
			legend(rows(1) stack size(small) ///
			symplacement(center) pos(6) ///
			label(1 "not at all meaningful") 	///
			label(2 "rather not meaningful")	///
			label(3 "neutral") 	///
			label(4 "rather meaningful") 	///
			label(5 "very meaningful")) 	
			graph save "dump/fig_S01c.gph", replace
			
		catplot incentive_money gender if gender < 3 & filter_1 == 1 & framing == 0, percent(gender) stack asyvars recast(bar) ///
			fxsize(25)									///
			title("Direct monetary" "compensation")	///
			ytitle("")									///
			l1title("", )								///
			ylabel(,)					///
			blabel(bar, pos(center) format(%2.1f) )	///
			var1opts(gap(*0) label(labsize(large))) ///
				bar(1, color(538t) fintensity(inten90)) ///
				bar(2, color(538t) fintensity(inten50)) ///
				bar(3, color(gray) fintensity(inten40)) ///
				bar(4, color(538o) fintensity(inten50)) ///
				bar(5, color(538o) fintensity(inten90)) ///
			legend(rows(1) stack size(small) ///
			symplacement(center) pos(6) ///
			label(1 "not at all meaningful") 	///
			label(2 "rather not meaningful")	///
			label(3 "neutral") 	///
			label(4 "rather meaningful") 	///
			label(5 "very meaningful")) 	
			graph save "dump/fig_S01d.gph", replace

			grc1leg 	"dump/fig_S01a.gph" ///
						"dump/fig_S01b.gph" ///
						"dump/fig_S01c.gph" ///
						"dump/fig_S01d.gph" , col(4) imargin (1 0 0 0) commonscheme title("Outcome I: How would you rate the following incentives" "to encourage people to get a COVID-19 vaccination?")

			graph save "dump/fig_S01A.gph", replace			
			
			catplot   gender money_incentive if gender < 3 & filter_1 == 1 & framing == 0, percent(gender) 				///
						asyvars recast(hbar) 												///
						legend(col(1) ring(0) position(5) bmargin(medium) size(vsmall))		///
						fysize(40) 															///
						ylabel(,labsize(vsmall))											///
						ytitle("percent" ,size(small))										///
						var2opts( label(labsize(vsmall)))								///
						title("Outcome II: If there were a financial incentive for COVID-19 vaccination," "how much money would you consider appropriate?") 
			graph save "dump/fig_S01B.gph", replace
			
			graph combine 	"dump/fig_S01A.gph" ///
							"dump/fig_S01B.gph", col(1) imargin (0 0 0 0) xsize(6) ysize(7)
							
				gr_edit			.AddTextBox added_text editor 50.83140785129991 -1.928044062301026
				gr_edit			.added_text_new = 1
				gr_edit			.added_text_rec = 1
				gr_edit			.added_text[1].style.editstyle  angle(default) size(medsmall) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
				gr_edit			.added_text[1].text = {}
				gr_edit			.added_text[1].text.Arrpush `"# SARS-CoV-2"'
				gr_edit			.added_text[1].text.Arrpush `"vaccination doses"'
				gr_edit			.added_text[1].style.editstyle size(tiny) editcopy
				gr_edit			.added_text[1].DragBy -1.087823750795866 1.813039584659828
				gr_edit			.added_text[1].DragBy -.1813039584659832 -.362607916931964
				gr_edit			.AddTextBox added_text editor 111.0709682260775 -.6897204236523002
				gr_edit			.added_text_new = 2
				gr_edit			.added_text_rec = 2
				gr_edit			.added_text[2].style.editstyle  angle(default) size(medlarge) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
				gr_edit			.added_text[2].text = {}
				gr_edit			.added_text[2].text.Arrpush a
				gr_edit			.AddTextBox added_text editor 38.3258352986633 -.8155770550146082
				gr_edit			.added_text_new = 3
				gr_edit			.added_text_rec = 3
				gr_edit			.added_text[3].style.editstyle  angle(default) size(medlarge) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
				gr_edit			.added_text[3].text = {}
				gr_edit			.added_text[3].text.Arrpush b
				gr_edit			.added_text[3].DragBy -.7551397881738351 .7551397881738504
			graph export "fig\fig_S01.emf", replace
			graph export "fig\fig_S01.pdf", replace


	
	********************************************************************************************************************************************************************	
	* Fig S02: Effects of frames – bivariate regression analyses (reference = control group; + 95% CI, full sample)
	********************************************************************************************************************************************************************

	coefplot m0_ologit_football_lw m0_ologit_sausage_lw m0_ologit_lottery_lw m0_ologit_money_lw , eform xline(1) keep(*framing*) drop(_cons) title(Dependent variable as ordinal scale) xtitle("Odds ratio from ordinal logit") legend(order(2 "Football ticket" 4 "Sausage" 6 "Lottery" 8 "Money" ) row(1)) 
	graph save "dump\m_S02a.gph", replace
	coefplot m0_ols_football_lw m0_ols_sausage_lw m0_ols_lottery_lw m0_ols_money_lw , xline(0) drop(_cons) title(Dependent variable used as metric scale) xtitle("b-coefficient from OLS") legend(order(2 "Football ticket" 4 "Sausage" 6 "Lottery" 8 "Money" ) row(1))
	graph save "dump\m_S02b.gph", replace
	coefplot m0_ologit_howmuchmoney_lw , pstyle(p8) eform xline(1) keep(*framing*) drop(_cons) title(Dependent variable as ordinal scale) xtitle("Odds ratio from ordinal logit") legend(order(2 "How much money" ) row(1))
	graph save "dump\m_S02c.gph", replace
	coefplot m0_ols_howmuchmoney_lw ,  pstyle(p8) xline(0) keep(*framing*) drop(_cons) title(Dependent variable rescaled to metric scale (in EUR)) xtitle("b-coefficient from OLS") legend(order(2 "How much money" ) row(1))
	graph save "dump\m_S02d.gph", replace

	grc1leg "dump\m_S02b.gph" "dump\m_S02a.gph", title ("Outcome I: Attitudes towards incentives") fysize(60)
	graph save "dump\m_S02_part1.gph", replace
	
	grc1leg "dump\m_S02d.gph" "dump\m_S02c.gph", title ("Outcome II: How much money should be paid as incentive to get vaccinated?") fysize(40)
	graph save "dump\m_S02_part2.gph", replace
	
	graph combine "dump\m_S02_part1.gph" "dump\m_S02_part2.gph", col(1) 														
	gr_edit			.plotregion1.graph1.plotregion1.graph1.AddTextBox added_text editor 39.411979885639 -1.831618735556752
	gr_edit			.plotregion1.graph1.plotregion1.graph1.added_text_new = 1
	gr_edit			.plotregion1.graph1.plotregion1.graph1.added_text_rec = 1
	gr_edit			.plotregion1.graph1.plotregion1.graph1.added_text[1].style.editstyle  angle(default) size(huge) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
	gr_edit			.plotregion1.graph1.plotregion1.graph1.added_text[1].text = {}
	gr_edit			.plotregion1.graph1.plotregion1.graph1.added_text[1].text.Arrpush a


	gr_edit			.plotregion1.graph1.plotregion1.graph2.AddTextBox added_text editor 39.73526079897225 .61866919997675
	gr_edit			.plotregion1.graph1.plotregion1.graph2.added_text_new = 1
	gr_edit			.plotregion1.graph1.plotregion1.graph2.added_text_rec = 1
	gr_edit			.plotregion1.graph1.plotregion1.graph2.added_text[1].style.editstyle  angle(default) size(huge) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
	gr_edit			.plotregion1.graph1.plotregion1.graph2.added_text[1].text = {}
	gr_edit			.plotregion1.graph1.plotregion1.graph2.added_text[1].text.Arrpush b


	gr_edit			.plotregion1.graph2.plotregion1.graph1.AddTextBox added_text editor 23.609400466656 -1.400577517779104
	gr_edit			.plotregion1.graph2.plotregion1.graph1.added_text_new = 1
	gr_edit			.plotregion1.graph2.plotregion1.graph1.added_text_rec = 1
	gr_edit			.plotregion1.graph2.plotregion1.graph1.added_text[1].style.editstyle  angle(default) size(huge) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
	gr_edit			.plotregion1.graph2.plotregion1.graph1.added_text[1].text = {}
	gr_edit			.plotregion1.graph2.plotregion1.graph1.added_text[1].text.Arrpush c

	gr_edit			.plotregion1.graph2.plotregion1.graph2.AddTextBox added_text editor 22.53179742221188 2.235073766642907
	gr_edit			.plotregion1.graph2.plotregion1.graph2.added_text_new = 1
	gr_edit			.plotregion1.graph2.plotregion1.graph2.added_text_rec = 1
	gr_edit			.plotregion1.graph2.plotregion1.graph2.added_text[1].style.editstyle  angle(default) size(huge) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
	gr_edit			.plotregion1.graph2.plotregion1.graph2.added_text[1].text = {}
	gr_edit			.plotregion1.graph2.plotregion1.graph2.added_text[1].text.Arrpush d


	gr_edit			.plotregion1.graph2.plotregion1.graph2.added_text[1].DragBy -.6465618266664661 -1.508644262221766


	gr_edit			.plotregion1.graph2.plotregion1.graph1.added_text[1].DragBy -1.724164871110595 -.1077603044444093


	graph export "fig\fig_S02.emf", replace
	graph export "fig\fig_S02.pdf", replace		
	
	
	
	
	********************************************************************************************************************************************************************	
	* Fig S03: Comparison OLS/Ologit
	********************************************************************************************************************************************************************	
	
	coefplot m1_ologit*, keep(*framing*) eform xline(1) drop(_cons) title(Dependent variable as ordinal scale) xtitle("Odds ratio from ordinal logit") legend(order(2 "Football" 4 "Sausage" 6 "Lottery" 8 "Money" ) row(1))
	graph save "dump\m_S03a.gph", replace
	coefplot m1_ols*, keep(*framing*) xline(0) drop(_cons) title(Dependent variable used as metric scale) xtitle("b-coefficient from OLS") legend(order(2 "Football" 4 "Sausage" 6 "Lottery" 8 "Money" ) row(1))
	graph save "dump\m_S03b.gph", replace
	coefplot m2_ologit_howmuchmoney , pstyle(p8) eform xline(1) keep(*framing*) drop(_cons) title(Dependent variable as ordinal scale) xtitle("Odds ratio from ordinal logit") legend(order(2 "How much money" ) row(1))
	graph save "dump\m_S03c.gph", replace
	coefplot m2_ols_howmuchmoney ,  pstyle(p8) xline(0) keep(*framing*) drop(_cons) title(Dependent variable rescaled to metric scale (in EUR)) xtitle("b-coefficient from OLS") legend(order(2 "How much money" ) row(1))
	graph save "dump\m_S03d.gph", replace

	grc1leg "dump\m_S03b.gph" "dump\m_S03a.gph", title ("Outcome I: Attitudes towards incentives") fysize(60)
	graph save "dump\m_S03_part1.gph", replace
	
	grc1leg "dump\m_S03d.gph" "dump\m_S03c.gph", title ("Outcome II: How much money should be paid as incentive to get vaccinated?") fysize(40)
	graph save "dump\m_S03_part2.gph", replace
	
	graph combine "dump\m_S03_part1.gph" "dump\m_S03_part2.gph", col(1) 	
	
		gr_edit			.plotregion1.graph1.plotregion1.graph1.AddTextBox added_text editor 39.411979885639 -1.831618735556752
		gr_edit			.plotregion1.graph1.plotregion1.graph1.added_text_new = 1
		gr_edit			.plotregion1.graph1.plotregion1.graph1.added_text_rec = 1
		gr_edit			.plotregion1.graph1.plotregion1.graph1.added_text[1].style.editstyle  angle(default) size(huge) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
		gr_edit			.plotregion1.graph1.plotregion1.graph1.added_text[1].text = {}
		gr_edit			.plotregion1.graph1.plotregion1.graph1.added_text[1].text.Arrpush a


		gr_edit			.plotregion1.graph1.plotregion1.graph2.AddTextBox added_text editor 39.73526079897225 .61866919997675
		gr_edit			.plotregion1.graph1.plotregion1.graph2.added_text_new = 1
		gr_edit			.plotregion1.graph1.plotregion1.graph2.added_text_rec = 1
		gr_edit			.plotregion1.graph1.plotregion1.graph2.added_text[1].style.editstyle  angle(default) size(huge) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
		gr_edit			.plotregion1.graph1.plotregion1.graph2.added_text[1].text = {}
		gr_edit			.plotregion1.graph1.plotregion1.graph2.added_text[1].text.Arrpush b


		gr_edit			.plotregion1.graph2.plotregion1.graph1.AddTextBox added_text editor 23.609400466656 -1.400577517779104
		gr_edit			.plotregion1.graph2.plotregion1.graph1.added_text_new = 1
		gr_edit			.plotregion1.graph2.plotregion1.graph1.added_text_rec = 1
		gr_edit			.plotregion1.graph2.plotregion1.graph1.added_text[1].style.editstyle  angle(default) size(huge) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
		gr_edit			.plotregion1.graph2.plotregion1.graph1.added_text[1].text = {}
		gr_edit			.plotregion1.graph2.plotregion1.graph1.added_text[1].text.Arrpush c

		gr_edit			.plotregion1.graph2.plotregion1.graph2.AddTextBox added_text editor 22.53179742221188 2.235073766642907
		gr_edit			.plotregion1.graph2.plotregion1.graph2.added_text_new = 1
		gr_edit			.plotregion1.graph2.plotregion1.graph2.added_text_rec = 1
		gr_edit			.plotregion1.graph2.plotregion1.graph2.added_text[1].style.editstyle  angle(default) size(huge) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
		gr_edit			.plotregion1.graph2.plotregion1.graph2.added_text[1].text = {}
		gr_edit			.plotregion1.graph2.plotregion1.graph2.added_text[1].text.Arrpush d


		gr_edit			.plotregion1.graph2.plotregion1.graph2.added_text[1].DragBy -.6465618266664661 -1.508644262221766


		gr_edit			.plotregion1.graph2.plotregion1.graph1.added_text[1].DragBy -1.724164871110595 -.1077603044444093
// 	gr_edit plotregion1.graph1.legend.AddTextBox added_text editor -.4743580015008895 2.856144699707229
// 	gr_edit plotregion1.graph1.legend.added_text_new = 1
// 	gr_edit plotregion1.graph1.legend.added_text_rec = 1
// 	gr_edit plotregion1.graph1.legend.added_text[1].style.editstyle  angle(default) size(medsmall) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
// 	gr_edit plotregion1.graph1.legend.added_text[1].style.editstyle size(vsmall) editcopy
// 	gr_edit plotregion1.graph1.legend.added_text[1].text = {}
// 	gr_edit plotregion1.graph1.legend.added_text[1].text.Arrpush n = 6,675
//
// 	gr_edit plotregion1.graph1.legend.AddTextBox added_text editor -.4743580015008895 18.43857995608003
// 	gr_edit plotregion1.graph1.legend.added_text_new = 2
// 	gr_edit plotregion1.graph1.legend.added_text_rec = 2
// 	gr_edit plotregion1.graph1.legend.added_text[2].style.editstyle  angle(default) size(medsmall) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
// 	gr_edit plotregion1.graph1.legend.added_text[2].style.editstyle size(vsmall) editcopy
// 	gr_edit plotregion1.graph1.legend.added_text[2].text = {}
// 	gr_edit plotregion1.graph1.legend.added_text[2].text.Arrpush n = 6,675
//
// 	gr_edit plotregion1.graph1.legend.AddTextBox added_text editor -.6455935537687297 34.53472186925628
// 	gr_edit plotregion1.graph1.legend.added_text_new = 3
// 	gr_edit plotregion1.graph1.legend.added_text_rec = 3
// 	gr_edit plotregion1.graph1.legend.added_text[3].style.editstyle  angle(default) size(medsmall) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
// 	gr_edit plotregion1.graph1.legend.added_text[3].style.editstyle size(vsmall) editcopy
// 	gr_edit plotregion1.graph1.legend.added_text[3].text = {}
// 	gr_edit plotregion1.graph1.legend.added_text[3].text.Arrpush n = 6,675
//
// 	gr_edit plotregion1.graph1.legend.AddTextBox added_text editor -.6455935537687297 48.91850825975428
// 	gr_edit plotregion1.graph1.legend.added_text_new = 4
// 	gr_edit plotregion1.graph1.legend.added_text_rec = 4
// 	gr_edit plotregion1.graph1.legend.added_text[4].style.editstyle  angle(default) size(medsmall) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
// 	gr_edit plotregion1.graph1.legend.added_text[4].text = {}
// 	gr_edit plotregion1.graph1.legend.added_text[4].text.Arrpush n = 6,675
// 	gr_edit plotregion1.graph1.legend.added_text[4].style.editstyle size(vsmall) editcopy
//
// 	gr_edit plotregion1.graph2.legend.AddTextBox added_text editor -.5129334319650113 2.951745333294005
// 	gr_edit plotregion1.graph2.legend.added_text_new = 1
// 	gr_edit plotregion1.graph2.legend.added_text_rec = 1
// 	gr_edit plotregion1.graph2.legend.added_text[1].style.editstyle  angle(default) size(medsmall) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
// 	gr_edit plotregion1.graph2.legend.added_text[1].style.editstyle size(vsmall) editcopy
// 	gr_edit plotregion1.graph2.legend.added_text[1].text = {}
// 	gr_edit plotregion1.graph2.legend.added_text[1].text.Arrpush n = 6,675

	graph export "fig\fig_S03.emf", replace
	graph export "fig\fig_S03.pdf", replace	

	
	********************************************************************************************************************************************************************	
	* Fig S04:  Predicted meaningfulness of vaccination incentives based on multiple OLS models with pooled frames
	********************************************************************************************************************************************************************		
	
	graph combine "dump/fig_football_ols_pooled2.gph" "dump/fig_sausage_ols_pooled2.gph" "dump/fig_lottery_ols_pooled2.gph" "dump/fig_money_ols_pooled2.gph", row(1) xsize(8) ysize(3) scale(1.5) title("Long COVID frames & Economic frames pooled") imargin(0 0 0 0 )
	gr_edit .plotregion1.graph1.xaxis1.major.num_rule_ticks = 0
	gr_edit .plotregion1.graph1.xaxis1.edit_tick 1 0 `"Control group"', custom tickset(major) editstyle(tickstyle(textstyle(color(reddish))) )
	gr_edit .plotregion1.graph2.xaxis1.major.num_rule_ticks = 0
	gr_edit .plotregion1.graph2.xaxis1.edit_tick 1 0 `"Control group"', custom tickset(major) editstyle(tickstyle(textstyle(color(reddish))) )
	gr_edit .plotregion1.graph3.xaxis1.major.num_rule_ticks = 0
	gr_edit .plotregion1.graph3.xaxis1.edit_tick 1 0 `"Control group"', custom tickset(major) editstyle(tickstyle(textstyle(color(reddish))) )
	gr_edit .plotregion1.graph4.xaxis1.major.num_rule_ticks = 0
	gr_edit .plotregion1.graph4.xaxis1.edit_tick 1 0 `"Control group"', custom tickset(major) editstyle(tickstyle(textstyle(color(reddish))) )
	gr_edit .plotregion1.graph1.plotregion1.AddMarker added_markers editor .0054740657037635 .2233700283787271
	gr_edit .plotregion1.graph1.plotregion1.added_markers_new = 1
	gr_edit .plotregion1.graph1.plotregion1.added_markers_rec = 1
	gr_edit .plotregion1.graph1.plotregion1.added_markers[1].style.editstyle  marker( symbol(square) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(reddish) size(medsmall) angle(stdarrow) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(zero) backangle(stdarrow) backsymangle(zero)) line( width(thin) color(black) pattern(solid) align(inside)) area( linestyle( width(thin) color(white) pattern(solid) align(inside)) shadestyle( color(white) intensity(inten100) fill(pattern10))) label( textstyle( horizontal(center) vertical(middle) angle(default) size(small) color(black) position() margin(zero) linestyle( width(thin) color(black) pattern(solid) align(inside))) position(6) textgap(half_tiny) horizontal(default) vertical(default)) dots( symbol(circle) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(black) size(vtiny) angle(horizontal) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(small) backangle(horizontal) backsymangle(zero)) connect(direct) connect_missings(yes) editcopy
	gr_edit .plotregion1.graph2.plotregion1.AddMarker added_markers editor -.009039258484574 .2064617743808659
	gr_edit .plotregion1.graph2.plotregion1.added_markers_new = 1
	gr_edit .plotregion1.graph2.plotregion1.added_markers_rec = 1
	gr_edit .plotregion1.graph2.plotregion1.added_markers[1].style.editstyle  marker( symbol(square) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(reddish) size(medsmall) angle(stdarrow) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(zero) backangle(stdarrow) backsymangle(zero)) line( width(thin) color(black) pattern(solid) align(inside)) area( linestyle( width(thin) color(white) pattern(solid) align(inside)) shadestyle( color(white) intensity(inten100) fill(pattern10))) label( textstyle( horizontal(center) vertical(middle) angle(default) size(small) color(black) position() margin(zero) linestyle( width(thin) color(black) pattern(solid) align(inside))) position(6) textgap(half_tiny) horizontal(default) vertical(default)) dots( symbol(circle) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(black) size(vtiny) angle(horizontal) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(small) backangle(horizontal) backsymangle(zero)) connect(direct) connect_missings(yes) editcopy
	gr_edit .plotregion1.graph3.plotregion1.AddMarker added_markers editor .0054740657037635 .2773912235717336
	gr_edit .plotregion1.graph3.plotregion1.added_markers_new = 1
	gr_edit .plotregion1.graph3.plotregion1.added_markers_rec = 1
	gr_edit .plotregion1.graph3.plotregion1.added_markers[1].style.editstyle  marker( symbol(square) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(reddish) size(medsmall) angle(stdarrow) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(zero) backangle(stdarrow) backsymangle(zero)) line( width(thin) color(black) pattern(solid) align(inside)) area( linestyle( width(thin) color(white) pattern(solid) align(inside)) shadestyle( color(white) intensity(inten100) fill(pattern10))) label( textstyle( horizontal(center) vertical(middle) angle(default) size(small) color(black) position() margin(zero) linestyle( width(thin) color(black) pattern(solid) align(inside))) position(6) textgap(half_tiny) horizontal(default) vertical(default)) dots( symbol(circle) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(black) size(vtiny) angle(horizontal) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(small) backangle(horizontal) backsymangle(zero)) connect(direct) connect_missings(yes) editcopy
	gr_edit .plotregion1.graph4.plotregion1.AddMarker added_markers editor -.009039258484574 .1905955197160464
	gr_edit .plotregion1.graph4.plotregion1.added_markers_new = 1
	gr_edit .plotregion1.graph4.plotregion1.added_markers_rec = 1
	gr_edit .plotregion1.graph4.plotregion1.added_markers[1].style.editstyle  marker( symbol(square) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(reddish) size(medsmall) angle(stdarrow) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(zero) backangle(stdarrow) backsymangle(zero)) line( width(thin) color(black) pattern(solid) align(inside)) area( linestyle( width(thin) color(white) pattern(solid) align(inside)) shadestyle( color(white) intensity(inten100) fill(pattern10))) label( textstyle( horizontal(center) vertical(middle) angle(default) size(small) color(black) position() margin(zero) linestyle( width(thin) color(black) pattern(solid) align(inside))) position(6) textgap(half_tiny) horizontal(default) vertical(default)) dots( symbol(circle) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(black) size(vtiny) angle(horizontal) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(small) backangle(horizontal) backsymangle(zero)) connect(direct) connect_missings(yes) editcopy

	graph save "dump/fig_pooled2", replace
	
	graph combine "dump/fig_football_ols_pooled3.gph" "dump/fig_sausage_ols_pooled3.gph" "dump/fig_lottery_ols_pooled3.gph" "dump/fig_money_ols_pooled3.gph", row(1) xsize(8) ysize(3) scale(1.5) title("All four frames pooled") imargin(0 0 0 0 )
	gr_edit .plotregion1.graph1.xaxis1.major.num_rule_ticks = 0
	gr_edit .plotregion1.graph1.xaxis1.edit_tick 1 0 `"Control group"', custom tickset(major) editstyle(tickstyle(textstyle(color(reddish))) )
	gr_edit .plotregion1.graph2.xaxis1.major.num_rule_ticks = 0
	gr_edit .plotregion1.graph2.xaxis1.edit_tick 1 0 `"Control group"', custom tickset(major) editstyle(tickstyle(textstyle(color(reddish))) )
	gr_edit .plotregion1.graph3.xaxis1.major.num_rule_ticks = 0
	gr_edit .plotregion1.graph3.xaxis1.edit_tick 1 0 `"Control group"', custom tickset(major) editstyle(tickstyle(textstyle(color(reddish))) )
	gr_edit .plotregion1.graph4.xaxis1.major.num_rule_ticks = 0
	gr_edit .plotregion1.graph4.xaxis1.edit_tick 1 0 `"Control group"', custom tickset(major) editstyle(tickstyle(textstyle(color(reddish))) )
	gr_edit .plotregion1.graph1.plotregion1.AddMarker added_markers editor .0054740657037635 .2233700283787271
	gr_edit .plotregion1.graph1.plotregion1.added_markers_new = 1
	gr_edit .plotregion1.graph1.plotregion1.added_markers_rec = 1
	gr_edit .plotregion1.graph1.plotregion1.added_markers[1].style.editstyle  marker( symbol(square) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(reddish) size(medsmall) angle(stdarrow) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(zero) backangle(stdarrow) backsymangle(zero)) line( width(thin) color(black) pattern(solid) align(inside)) area( linestyle( width(thin) color(white) pattern(solid) align(inside)) shadestyle( color(white) intensity(inten100) fill(pattern10))) label( textstyle( horizontal(center) vertical(middle) angle(default) size(small) color(black) position() margin(zero) linestyle( width(thin) color(black) pattern(solid) align(inside))) position(6) textgap(half_tiny) horizontal(default) vertical(default)) dots( symbol(circle) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(black) size(vtiny) angle(horizontal) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(small) backangle(horizontal) backsymangle(zero)) connect(direct) connect_missings(yes) editcopy
	gr_edit .plotregion1.graph2.plotregion1.AddMarker added_markers editor -.009039258484574 .2064617743808659
	gr_edit .plotregion1.graph2.plotregion1.added_markers_new = 1
	gr_edit .plotregion1.graph2.plotregion1.added_markers_rec = 1
	gr_edit .plotregion1.graph2.plotregion1.added_markers[1].style.editstyle  marker( symbol(square) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(reddish) size(medsmall) angle(stdarrow) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(zero) backangle(stdarrow) backsymangle(zero)) line( width(thin) color(black) pattern(solid) align(inside)) area( linestyle( width(thin) color(white) pattern(solid) align(inside)) shadestyle( color(white) intensity(inten100) fill(pattern10))) label( textstyle( horizontal(center) vertical(middle) angle(default) size(small) color(black) position() margin(zero) linestyle( width(thin) color(black) pattern(solid) align(inside))) position(6) textgap(half_tiny) horizontal(default) vertical(default)) dots( symbol(circle) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(black) size(vtiny) angle(horizontal) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(small) backangle(horizontal) backsymangle(zero)) connect(direct) connect_missings(yes) editcopy
	gr_edit .plotregion1.graph3.plotregion1.AddMarker added_markers editor .0054740657037635 .2773912235717336
	gr_edit .plotregion1.graph3.plotregion1.added_markers_new = 1
	gr_edit .plotregion1.graph3.plotregion1.added_markers_rec = 1
	gr_edit .plotregion1.graph3.plotregion1.added_markers[1].style.editstyle  marker( symbol(square) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(reddish) size(medsmall) angle(stdarrow) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(zero) backangle(stdarrow) backsymangle(zero)) line( width(thin) color(black) pattern(solid) align(inside)) area( linestyle( width(thin) color(white) pattern(solid) align(inside)) shadestyle( color(white) intensity(inten100) fill(pattern10))) label( textstyle( horizontal(center) vertical(middle) angle(default) size(small) color(black) position() margin(zero) linestyle( width(thin) color(black) pattern(solid) align(inside))) position(6) textgap(half_tiny) horizontal(default) vertical(default)) dots( symbol(circle) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(black) size(vtiny) angle(horizontal) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(small) backangle(horizontal) backsymangle(zero)) connect(direct) connect_missings(yes) editcopy
	gr_edit .plotregion1.graph4.plotregion1.AddMarker added_markers editor -.009039258484574 .1905955197160464
	gr_edit .plotregion1.graph4.plotregion1.added_markers_new = 1
	gr_edit .plotregion1.graph4.plotregion1.added_markers_rec = 1
	gr_edit .plotregion1.graph4.plotregion1.added_markers[1].style.editstyle  marker( symbol(square) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(reddish) size(medsmall) angle(stdarrow) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(zero) backangle(stdarrow) backsymangle(zero)) line( width(thin) color(black) pattern(solid) align(inside)) area( linestyle( width(thin) color(white) pattern(solid) align(inside)) shadestyle( color(white) intensity(inten100) fill(pattern10))) label( textstyle( horizontal(center) vertical(middle) angle(default) size(small) color(black) position() margin(zero) linestyle( width(thin) color(black) pattern(solid) align(inside))) position(6) textgap(half_tiny) horizontal(default) vertical(default)) dots( symbol(circle) linestyle( width(thin) color(black) pattern(solid) align(inside)) fillcolor(black) size(vtiny) angle(horizontal) symangle(zero) backsymbol(none) backline( width(thin) color(black) pattern(solid) align(inside)) backcolor(black) backsize(small) backangle(horizontal) backsymangle(zero)) connect(direct) connect_missings(yes) editcopy

	graph save "dump/fig_pooled3", replace	
	
	graph combine "dump/fig_pooled2" "dump/fig_pooled3", col(1) xsize(6) ysize(6) imargin(0 0 0 0 ) scale(0.8)
	graph export "fig\fig_S04.emf", replace			
	graph export "fig\fig_S04.pdf", replace	
	
	
	
	********************************************************************************************************************************************************************	
	* Fig S05: Descriptive statistics for concrete drivers of vaccination. Means by sex, age_group and party
	********************************************************************************************************************************************************************	
	
	use "dump\PPD_July_2022_labeled.dta", replace
	graph bar (mean) vacc_reason_1 if gender <3 & filter_1 == 1,  over(age_group_r) asyvars by(gender, row(1) title("Protecting self") note("")) ytitle(Importance)  caption("") legend(rows(1) )  ysc(r(0 1)) 
	graph save "dump\S05_1a.gph", replace
	graph bar (mean) vacc_reason_2 if gender <3 & filter_1 == 1,  over(age_group_r) asyvars by(gender, row(1) title("Protecting others") note("")) ytitle(Importance)  caption("") legend(rows(1) ) ysc(r(0 1))
	graph save "dump\S05_1b.gph", replace
	graph bar (mean) vacc_reason_3 if gender <3 & filter_1 == 1,  over(age_group_r) asyvars by(gender, row(1) title("Participation in public events") note("")) ytitle(Importance)  caption("") legend(rows(1) ) ysc(r(0 1))
	graph save "dump\S05_1c.gph", replace
	graph bar (mean) vacc_reason_4 if gender <3 & filter_1 == 1,  over(age_group_r) asyvars by(gender, row(1) title("Vocational mandates") note("")) ytitle(Importance)  caption("") legend(rows(1) ) ysc(r(0 1))
	graph save "dump\S05_1d.gph", replace
	graph bar (mean) vacc_reason_5 if gender <3 & filter_1 == 1,  over(age_group_r) asyvars by(gender, row(1) title("Medical advice") note("")) ytitle(Importance)  caption("") legend(rows(1) ) ysc(r(0 1))
	graph save "dump\S05_1e.gph", replace
	graph bar (mean) vacc_reason_6 if gender <3 & filter_1 == 1,  over(age_group_r) asyvars by(gender, row(1) title("Peer pressure") note("")) ytitle(Importance)  caption("") legend(rows(1) ) ysc(r(0 1))
	graph save "dump\S05_1f.gph", replace
	
	grc1leg2 "dump\S05_1a.gph" "dump\S05_1b.gph" "dump\S05_1c.gph" "dump\S05_1d.gph" "dump\S05_1e.gph" "dump\S05_1f.gph", xcommon
	graph save "dump\S05_part1.gph", replace
	
	graph bar (mean) vacc_reason_1 if  filter_1 == 1,  over(sonntagsfrage) asyvars bar(1, lcolor(black) fcolor(black)) bar(2, lcolor(red) fcolor(red)) bar(3, lcolor(green) fcolor(green)) bar(4, lcolor(yellow) fcolor(yellow)) bar(5, lcolor(purple) fcolor(purple)) bar(6, lcolor(eltblue) fcolor(eltblue)) bar(7, lcolor(gray) fcolor(gray)) title("Protecting self")  ytitle(Importance)   legend(rows(1))
	graph save "dump\S05_2a.gph", replace
	graph bar (mean) vacc_reason_2 if  filter_1 == 1,  over(sonntagsfrage) asyvars bar(1, lcolor(black) fcolor(black)) bar(2, lcolor(red) fcolor(red)) bar(3, lcolor(green) fcolor(green)) bar(4, lcolor(yellow) fcolor(yellow)) bar(5, lcolor(purple) fcolor(purple)) bar(6, lcolor(eltblue) fcolor(eltblue)) bar(7, lcolor(gray) fcolor(gray)) title("Protecting others")  ytitle(Importance)   legend(rows(1))
	graph save "dump\S05_2b.gph", replace
	graph bar (mean) vacc_reason_3 if  filter_1 == 1,  over(sonntagsfrage) asyvars bar(1, lcolor(black) fcolor(black)) bar(2, lcolor(red) fcolor(red)) bar(3, lcolor(green) fcolor(green)) bar(4, lcolor(yellow) fcolor(yellow)) bar(5, lcolor(purple) fcolor(purple)) bar(6, lcolor(eltblue) fcolor(eltblue)) bar(7, lcolor(gray) fcolor(gray)) title("Participation in public events")  ytitle(Importance)   legend(rows(1))
	graph save "dump\S05_2c.gph", replace
	graph bar (mean) vacc_reason_4 if  filter_1 == 1,  over(sonntagsfrage) asyvars bar(1, lcolor(black) fcolor(black)) bar(2, lcolor(red) fcolor(red)) bar(3, lcolor(green) fcolor(green)) bar(4, lcolor(yellow) fcolor(yellow)) bar(5, lcolor(purple) fcolor(purple)) bar(6, lcolor(eltblue) fcolor(eltblue)) bar(7, lcolor(gray) fcolor(gray)) title("Vocational mandates")  ytitle(Importance)   legend(rows(1))
	graph save "dump\S05_2d.gph", replace
	graph bar (mean) vacc_reason_5 if  filter_1 == 1,  over(sonntagsfrage) asyvars bar(1, lcolor(black) fcolor(black)) bar(2, lcolor(red) fcolor(red)) bar(3, lcolor(green) fcolor(green)) bar(4, lcolor(yellow) fcolor(yellow)) bar(5, lcolor(purple) fcolor(purple)) bar(6, lcolor(eltblue) fcolor(eltblue)) bar(7, lcolor(gray) fcolor(gray)) title("Medical advice")  ytitle(Importance)   legend(rows(1))
	graph save "dump\S05_2e.gph", replace
	graph bar (mean) vacc_reason_6 if  filter_1 == 1,  over(sonntagsfrage) asyvars bar(1, lcolor(black) fcolor(black)) bar(2, lcolor(red) fcolor(red)) bar(3, lcolor(green) fcolor(green)) bar(4, lcolor(yellow) fcolor(yellow)) bar(5, lcolor(purple) fcolor(purple)) bar(6, lcolor(eltblue) fcolor(eltblue)) bar(7, lcolor(gray) fcolor(gray)) title("Peer pressure")  ytitle(Importance)   legend(rows(1))
	graph save "dump\S05_2f.gph", replace
	
	grc1leg2 "dump\S05_2a.gph" "dump\S05_2b.gph" "dump\S05_2c.gph" "dump\S05_2d.gph" "dump\S05_2e.gph" "dump\S05_2f.gph", ycommon
	graph save "dump\S05_part2.gph", replace
	
	graph combine "dump\S05_part1.gph" "dump\S05_part2.gph", rows(2) xsize(5) ysize(6.5) 
	gr_edit .plotregion1.graph1.AddTextBox added_text editor 57.02124153642194 -.0220806191107501
	gr_edit .plotregion1.graph1.added_text_new = 1
	gr_edit .plotregion1.graph1.added_text_rec = 1
	gr_edit .plotregion1.graph1.added_text[1].style.editstyle  angle(default) size(medsmall) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
	gr_edit .plotregion1.graph1.added_text[1].style.editstyle size(medlarge) editcopy
	gr_edit .plotregion1.graph1.added_text[1].text = {}
	gr_edit .plotregion1.graph1.added_text[1].text.Arrpush `"a"'
	gr_edit .plotregion1.graph2.AddTextBox added_text editor 56.36611641475782 -.3165359561022312
	gr_edit .plotregion1.graph2.added_text_new = 1
	gr_edit .plotregion1.graph2.added_text_rec = 1
	gr_edit .plotregion1.graph2.added_text[1].style.editstyle  angle(default) size(medlarge) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
	gr_edit .plotregion1.graph2.added_text[1].text = {}
	gr_edit .plotregion1.graph2.added_text[1].text.Arrpush `"b"'
	gr_edit .plotregion1.graph2.added_text[1].DragBy 1.0305936794702 .5889106739829616
	graph export "fig\fig_S05.emf", replace
	graph export "fig\fig_S05.pdf", replace 
	
	
	

	
	********************************************************************************************************************************************************************	
	* Fig S06: Distribution of main demographic variables in the sample compared to reality
	********************************************************************************************************************************************************************
 	
	* Read data 
	use "dump\PPD_July_2022_labeled.dta", replace
	* Gender (data from Destatis 12411-0013:  Bevölkerung: Bundesländer, Stichtag, Geschlecht Fortschreibung des Bevölkerungsstandes)
	gen		str gender2 	= "male" in 1
		replace gender2 	= "female" in 2
		replace gender2 	= "non-binary" in 3
		
	gen gender_real			= 48.9250541 in 1
		replace gender_real	= 51.0749459 in 2
		replace gender_real	= 0 in 3

	estpost tab gender if filter_1 == 1
	mat gender_survey = e(pct)'	
	svmat gender_survey	
	rename gender_survey1 gender_survey
	replace gender_survey = . if gender_real == .

	label variable gender_real "Official statistics, Dec 2021"
	label variable gender_survey "Survey, Jun/Jul 2022"	
	graph bar (asis) gender_real gender_survey, over(gender2) bar(1, color(gs10) fint(80)) bar(2, color(538m))  legend(col(2) pos(6)) ///
	///note("N=6,675", ring(0) pos(1)) ///
	ytitle("Percent") title("Gender") 
	graph save "dump/gender_survey.gph", replace 

	* Age (data from Destatis 12411-0013:  Bevölkerung: Bundesländer, Stichtag, Geschlecht; Fortschreibung des Bevölkerungsstandes)
	gen age2 = _n +17
	replace age2 = . if age2 > 90

	gen age_real = .
	replace age_real = 0.0111418039055486 if age2 == 18
	replace age_real = 0.0114986529869714 if age2 == 19
	replace age_real = 0.0119848447250272 if age2 == 20
	replace age_real = 0.0127466878196854 if age2 == 21
	replace age_real = 0.0130162994378358 if age2 == 22
	replace age_real = 0.0134075274600889 if age2 == 23
	replace age_real = 0.014035285478184 if age2 == 24
	replace age_real = 0.0139619581523964 if age2 == 25
	replace age_real = 0.0137480735720865 if age2 == 26
	replace age_real = 0.0139275503822657 if age2 == 27
	replace age_real = 0.0143372868154311 if age2 == 28
	replace age_real = 0.0145521083480068 if age2 == 29
	replace age_real = 0.0150259899747549 if age2 == 30
	replace age_real = 0.0163484332320248 if age2 == 31
	replace age_real = 0.0161257413004162 if age2 == 32
	replace age_real = 0.0164583449401298 if age2 == 33
	replace age_real = 0.0160947498023932 if age2 == 34
	replace age_real = 0.0157837825526947 if age2 == 35
	replace age_real = 0.0152987439866584 if age2 == 36
	replace age_real = 0.0151791023896391 if age2 == 37
	replace age_real = 0.0151818411731853 if age2 == 38
	replace age_real = 0.0154614853879051 if age2 == 39
	replace age_real = 0.0153166181529601 if age2 == 40
	replace age_real = 0.0153750983889971 if age2 == 41
	replace age_real = 0.0146256086495974 if age2 == 42
	replace age_real = 0.0144086393341354 if age2 == 43
	replace age_real = 0.014254330503281 if age2 == 44
	replace age_real = 0.0139816341499785 if age2 == 45
	replace age_real = 0.0136046045582151 if age2 == 46
	replace age_real = 0.0137549493602526 if age2 == 47
	replace age_real = 0.0138373146717427 if age2 == 48
	replace age_real = 0.0149401795618566 if age2 == 49
	replace age_real = 0.0163708624278033 if age2 == 50
	replace age_real = 0.0169336680319022 if age2 == 51
	replace age_real = 0.0181318570040749 if age2 == 52
	replace age_real = 0.0189329656059959 if age2 == 53
	replace age_real = 0.0193634591355116 if age2 == 54
	replace age_real = 0.0198132106377524 if age2 == 55
	replace age_real = 0.0198048357259611 if age2 == 56
	replace age_real = 0.0200999901043426 if age2 == 57
	replace age_real = 0.0198443174529774 if age2 == 58
	replace age_real = 0.0191318301207522 if age2 == 59
	replace age_real = 0.0187465985930004 if age2 == 60
	replace age_real = 0.0180294841580471 if age2 == 61
	replace age_real = 0.0174170777424611 if age2 == 62
	replace age_real = 0.0163504945270096 if age2 == 63
	replace age_real = 0.015858392782354 if age2 == 64
	replace age_real = 0.0153056053601742 if age2 == 65
	replace age_real = 0.0147086082056982 if age2 == 66
	replace age_real = 0.0142985834795279 if age2 == 67
	replace age_real = 0.0137215362009887 if age2 == 68
	replace age_real = 0.0135766401367431 if age2 == 69
	replace age_real = 0.0131612387460321 if age2 == 70
	replace age_real = 0.0130163859257373 if age2 == 71
	replace age_real = 0.0124278069270034 if age2 == 72
	replace age_real = 0.0111877001519232 if age2 == 73
	replace age_real = 0.0103634127924111 if age2 == 74
	replace age_real = 0.00885357908197849 if age2 == 75
	replace age_real = 0.00764052860540493 if age2 == 76
	replace age_real = 0.0100166827954591 if age2 == 77
	replace age_real = 0.00996342066281012 if age2 == 78
	replace age_real = 0.00941658649118079 if age2 == 79
	replace age_real = 0.0110886859194021 if age2 == 80
	replace age_real = 0.0112094662737906 if age2 == 81
	replace age_real = 0.0105624791122709 if age2 == 82
	replace age_real = 0.00929123669266517 if age2 == 83
	replace age_real = 0.00810064424117065 if age2 == 84
	replace age_real = 0.00721644959524743 if age2 == 85
	replace age_real = 0.00631370329446111 if age2 == 86
	replace age_real = 0.00523884606976993 if age2 == 87
	replace age_real = 0.00369457575990613 if age2 == 88
	replace age_real = 0.00318010247807297 if age2 == 89
	replace age_real = 0.0122011797958785 if age2 == 90

	gen age_survey = age
	label variable age_real "Official statistics, Dec 2021"
	label variable age_survey "Survey, Jun/Jul 2022"

	twoway bar age_real age2,  lcolor(gs10) fcolor(gs10) || hist age_survey if filter_1 == 1, width(1)  fcolor(none) lcolor(538m) || pcarrowi 0.012 94 0.0118 91  , ///
	lcolor(gs6) mcolor(gs6) msize(vsmall) legend(order(1 2 "Survey, June 2022" ) col(2) pos(6)) ///
	///note("N=6,675", ring(0) pos(1)) ///
	ytitle("Density") xtitle("") text(0.0121 96  "90+", color(gs6) size(small)) title("Age")
	graph save "dump/age_survey.gph", replace

	* education (data from Destatis 12211-9012 Allgemeiner Schulabschluss & 12211-9015 Beruflicher Schulabschluss)
	gen			edu2 	= 1 in 1
		replace edu2 	= 2 in 2
		replace edu2 	= 3 in 3
		
	label values edu2 edu_3 
		
	gen edu_real			= 32.68595569 in 1 		// ohne allg. Schulabschluss + Haupt- (Volks-)schulabschluss
		replace edu_real	= 30.23489885 in 2		// polytechnische Oberschule + Realschulabschluss + ohne Angabe
		replace edu_real	= 37.07914546 in 3		// Abitur/Fachabitur

	estpost tab edu_3
	mat edu_survey = e(pct)'	
	svmat edu_survey	
	rename edu_survey1 edu_survey	
	replace edu_survey = . if edu_real == .

	label variable edu_real "Official statistics, Dec 2019"
	label variable edu_survey "Survey, Jun/Jul 2022"	
	*splitvallabels edu2 , length(20)
	graph bar (asis) edu_real edu_survey, over(edu2) bar(1, color(gs10) fint(80)) bar(2, color(538m)) legend(col(2) pos(6)) ///
	///note("N=6,675", ring(0) pos(1)) ///
	ytitle("Percent") title("Education")
	graph save "dump/edu_survey.gph", replace
	 
	* federal state (data from Destatis 12411-0013:  Bevölkerung: Bundesländer, Stichtag, Geschlecht)
	gen		bundesland2 	= _n
		replace bundesland2 = . if bundesland2 > 16
		label values bundesland2 bundesland

	gen bundesland_real = .
	replace bundesland_real = 13.29371947 if bundesland2 == 1
	replace bundesland_real = 15.82356122 if bundesland2 == 2
	replace bundesland_real = 4.410665313 if bundesland2 == 3
	replace bundesland_real = 3.069779952 if bundesland2 == 4
	replace bundesland_real = 0.811429491 if bundesland2 == 5
	replace bundesland_real = 2.217066903 if bundesland2 == 6
	replace bundesland_real = 7.533469557 if bundesland2 == 7
	replace bundesland_real = 1.965747476 if bundesland2 == 8
	replace bundesland_real = 9.623673122 if bundesland2 == 9
	replace bundesland_real = 21.45456218 if bundesland2 == 10
	replace bundesland_real = 4.942643746 if bundesland2 == 11
	replace bundesland_real = 1.203162027 if bundesland2 == 12
	replace bundesland_real = 4.892636442 if bundesland2 == 13
	replace bundesland_real = 2.658832689 if bundesland2 == 14
	replace bundesland_real = 3.526124716 if bundesland2 == 15
	replace bundesland_real = 2.572925698 if bundesland2 == 16
		
	estpost tab bundesland
	mat bundesland_survey = e(pct)'	
	svmat bundesland_survey
	rename bundesland_survey1 bundesland_survey	
	replace bundesland_survey = . if bundesland_real == .

	label variable bundesland_real "Official statistics, Dec 2021"
	label variable bundesland_survey "Survey, Jun/Jul 2022"	
	graph hbar (asis) bundesland_real bundesland_survey, over(bundesland2) bar(1, color(gs10) fint(80)) bar(2, color(538m)) legend(col(2) pos(6)) ///
	///note("N=6,675", ring(0) pos(1)) ///
	ytitle("Percent") title("State")
	graph save "dump/bundesland_survey.gph", replace

	* party (mean of the following four opinion polls: 
			* 28.June-30.June 2022	Forschungsgruppe Wahlen
			* 28.June-02.July 2022 FORSA
			* 04.July-06.July 2022 Infratest Dimap
			* 01.July-05.July 2022	YouGov
			* source: https://www.wahlrecht.de/umfragen/
			
	gen party = sonntagsfrage
	replace party = . if sonntagsfrage > 7

	gen party2 = _n
	replace party2 = . if party2 > 7
	label values party2 sonntagsfrage

	gen party_real = .
	replace party_real = 27 if party2 == 1
	replace party_real = 19.75 if party2 == 2
	replace party_real = 23 if party2 == 3
	replace party_real = 7.25 if party2 == 4
	replace party_real = 4.5 if party2 == 5
	replace party_real = 10.5 if party2 == 6
	replace party_real = 8 if party2 == 7

	estpost tab party
	mat party_survey = e(pct)'	
	svmat party_survey	
	rename party_survey1 party_survey	
	replace party_survey = . if party_real == .

	label variable party_real "Opinion polls 28.Jun - 05.Jul 2022"
	label variable party_survey "Survey, Jun/Jul 2022"	
	graph bar (asis) party_real party_survey, over(party2) bar(1, color(gs10) fint(80)) bar(2, color(538m)) legend(col(2) pos(6)) ///
	///note("N=6,675", ring(0) pos(1)) ///
	ytitle("Percent") title("Voting intention")
	graph save "dump/party_survey.gph", replace

		graph combine "dump/gender_survey.gph" "dump/age_survey.gph" "dump/edu_survey.gph" "dump/party_survey.gph" ,  rows(2) xsize(5) ysize(5)
		graph save "dump/distribution_survey_a.gph", replace

		graph combine "dump/distribution_survey_a.gph" "dump/bundesland_survey.gph", rows(2) xsize(6) ysize(8)
			gr_edit	.plotregion1.graph2.gmetric_mult = .8
			gr_edit .plotregion1.graph1.plotregion1.graph4.legend.plotregion1.key[1].DragBy 0 -3.017241379310334
			// key[1] reposition
			gr_edit .plotregion1.graph1.plotregion1.graph4.legend.plotregion1.label[1].DragBy 0 -3.879310344827581
			// label[1] reposition
			gr_edit .gmetric_mult = .87
	graph export "fig\fig_S06.emf", replace
	graph export "fig\fig_S06.pdf", replace
	
	********************************************************************************************************************************************************************	
	* Fig S07: Money deemed acceptable as incentive for COVID-19 vaccination (pooled frames)
	********************************************************************************************************************************************************************		
	
	use "dump\PPD_July_2022_labeled.dta", replace
	reg money_incentive_metric 													/// Analysis II DV 
			i.framing_pool_types															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 

	margins , at(framing=(0(1)2)  )
	marginsplot, recast(scatter) title("") ytitle("") xtitle("Money to be paid as incentive for COVID-19 vaccination (in EUR)") horiz xsize(5) ysize(1.2) scale(2) level(95)
	graph save "dump\fig_07a.gph", replace
	
		reg money_incentive_metric 													/// Analysis II DV 
			i.framing_pool_all															/// frames
			i.vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
			threat_covid safety_perception 										/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  i.religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN i.sonntagsfrage spaltungsindex					/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
			b10.bundesland i.gender i.age_group_r i.edu_3 i.hh_inc  health_01   /// controls: sociodemographic    
			if filter_1 == 1 

	margins , at(framing=(0(1)1)  )
	marginsplot, recast(scatter) title("") ytitle("") xtitle("Money to be paid as incentive for COVID-19 vaccination (in EUR)") horiz xsize(5) ysize(1.2) scale(2) level(95)
	graph save "dump\fig_07b.gph", replace
	
	graph combine "dump\fig_07a.gph" "dump\fig_07b.gph", col(1) imargin( 0 0 0 0 ) scale(0.8)
	gr_edit .AddTextBox added_text editor 93.88556686017222 -.2980163721920903
	gr_edit .added_text_new = 1
	gr_edit .added_text_rec = 1
	gr_edit .added_text[1].style.editstyle  angle(default) size(half_tiny) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
	gr_edit .added_text[1].style.editstyle size(huge) editcopy
	gr_edit .added_text[1].text = {}
	gr_edit .added_text[1].text.Arrpush `"a)"'
// editor text[1] edits

	gr_edit .AddTextBox added_text editor 51.05249940344796 .2642850231186449
	gr_edit .added_text_new = 2
	gr_edit .added_text_rec = 2
	gr_edit .added_text[2].style.editstyle  angle(default) size(half_tiny) color(black) horizontal(left) vertical(middle) margin(zero) linegap(zero) drawbox(no) boxmargin(zero) fillcolor(bluishgray) linestyle( width(vthin) color(black) pattern(solid) align(inside)) box_alignment(east) editcopy
	gr_edit .added_text[2].style.editstyle size(huge) editcopy
	gr_edit .added_text[2].text = {}
	gr_edit .added_text[2].text.Arrpush `"b)"'
// editor text[2] edits
	graph export "fig\fig_S07.emf", replace
	graph export "fig\fig_S07.pdf", replace
	
******************************************************************
**																**
**		8)	Tables Online Annex									**
**																**
******************************************************************	


asdoc codebook, replace save(dump\codebook_stata.rtf) title(Codebook Politikpanel Deutschland Umfrage Juli 2022) 
* Word File mit Tabellen/Statistiken aus unicode umkodieren
unicode convertfile  dump\codebook_stata.rtf  dump\codebook_stata_unicode.rtf, dstencoding(Windows-1252) replace	

********************************************************************************************************************************************************************	
* Table S1: Means of outcome variables of control group and the pooled framing groups (+ results of mean differences tests)
********************************************************************************************************************************************************************
* 1) Dunnetts test control vs Long COVID/Economic frames	
		* Read data 
	use "dump\PPD_July_2022_labeled.dta", replace
	drop if filter_1 == 0
	bysort framing_pool_types: sum(incentive_football)
	pwmean incentive_football , over(framing_pool_types)  mcompare(dunnett) effects
	** --> no significant differences between experimental groups and control group
	
	bysort framing_pool_types: sum(incentive_sausage)
	pwmean incentive_sausage , over(framing_pool_types)  mcompare(dunnett) effects
	** --> no significant differences between experimental groups and control group
	
	bysort framing_pool_types: sum(incentive_lottery)	
	pwmean incentive_lottery , over(framing_pool_types)  mcompare(dunnett) effects
	** --> Long Covid Frames: p < 0.05; Economic Frames: p < 0.05
	
	bysort framing_pool_types: sum(incentive_money)	
	pwmean incentive_money , over(framing_pool_types)  mcompare(dunnett) effects
	** --> Long Covid Frames: p < 0.05; Economic Frames: p < 0.1

	bysort framing_pool_types: sum(money_incentive_metric)	
	pwmean money_incentive_metric , over(framing_pool_types)  mcompare(dunnett) effects
	** --> no significance

	*2) t-test control vs. pooled all frames 
	* Levene test for equal variances
	sdtest incentive_football, by(framing_pool_all)
	* --> equal variances
	sdtest incentive_sausage, by(framing_pool_all)
	* --> equal variances
	sdtest incentive_lottery, by(framing_pool_all)
	* --> equal variances
	sdtest incentive_money, by(framing_pool_all)
	* --> unequal variances
	
	ttest incentive_football, by(framing_pool_all)
	* --> no significant mean difference 
	ttest incentive_sausage, by(framing_pool_all)
	* --> no significant mean difference 
	ttest incentive_lottery, by(framing_pool_all)	
	* --> significant difference at p < 0.01
	ttest incentive_money, by(framing_pool_all) unequal
	* --> significant difference at p < 0.01
	
	bysort framing_pool_all: sum(money_incentive_metric)	
	sdtest money_incentive_metric, by(framing_pool_all)
	* --> unequal variances
	ttest money_incentive_metric , by(framing_pool_all)  unequal
	** --> --> significant difference at p < 0.1

	********************************************************************************************************************************************************************	
	* Table S2: Outcome I - Main models, OLS regression (b-coefficients)
	********************************************************************************************************************************************************************
	
	* see above
	
	********************************************************************************************************************************************************************	
	* Table S3: Outcome I - Cross check models A: concrete drivers for vaccination decision, OLS regression (b-coefficients)
	********************************************************************************************************************************************************************
	
	* see above
	
	********************************************************************************************************************************************************************	
	* Table S4: Outcome I - Cross check models B: self perceived knowledge about COVID-19 vaccination, OLS regression (b-coefficients)
	********************************************************************************************************************************************************************
	
	* see above
	
	********************************************************************************************************************************************************************	
	* Table S5: Outcome I - Cross check models C: private health insurance, OLS regression (b-coefficients)
	********************************************************************************************************************************************************************
	
	* see above
	
	********************************************************************************************************************************************************************	
	* Table S6: Outcome II - Main model & cross check models, OLS regression (b-coefficients)
	********************************************************************************************************************************************************************

	* see above


	********************************************************************************************************************************************************************	
	* Table S7: Survey Questions
	********************************************************************************************************************************************************************	


	********************************************************************************************************************************************************************	
	* Table S8: Missing values – item non-response
	********************************************************************************************************************************************************************	

		asdoc mdesc incentive_* 													/// Analysis I DV 
			money_incentive	money_incentive_metric									/// Analysis II DV
			framing 														 		///	main IV
			vacc_04 vacc_decision_* 												/// controls: COVID-19 vaccination 
			threat_covid safety_perception 											/// controls: perception COVID-19
			school_system_3  medical_system_2 medical_system_1  religion 			/// controls: attitudes towards esoteric and medical belief systems & religion
			left_right GAL_TAN sonntagsfrage spaltungsindex							/// controls: political attitudes  
			solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  			/// controls: psychosocial 
			bundesland gender age age_group_r edu_3 hh_inc  health_01   				/// controls: sociodemographic    
			informed_vaccination private_health_insurance vacc_reason* ,	label	/// controls: additional cross checks
			replace save(dump\table_S07_not_unicode.rtf) title(Table S07: Missing values – item non-response) 
			unicode convertfile dump\table_S07_not_unicode.rtf  tables\table_S07.rtf, dstencoding(Windows-1252) replace

	********************************************************************************************************************************************************************	
	* Table S9: Descriptive statistics for the original dataset and the final dataset used in the regression models (reduced via listwise deletion)
	********************************************************************************************************************************************************************	

	clear all
			
	* Read data 
	use "dump\PPD_July_2022_labeled.dta", replace	
	
	eststo: estpost summarize incentive_* 									/// Analysis I DV 
		money_incentive money_incentive_metric								/// Analysis II DV
		framing 														 	///	main IV
		vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
		threat_covid safety_perception 										/// controls: perception COVID-19
		school_system_3  medical_system_2 medical_system_1  religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
		left_right GAL_TAN sonntagsfrage spaltungsindex						/// controls: political attitudes  
		solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
		bundesland gender age_group_r edu_3 hh_inc  health_01   			/// controls: sociodemographic    
		informed_vaccination private_health_insurance  vacc_reason_*	
	
	eststo: estpost summarize incentive_* 									/// Analysis I DV 
		money_incentive money_incentive_metric								/// Analysis II DV
		framing 														 	///	main IV
		vacc_04 vacc_decision_* 											/// controls: COVID-19 vaccination 
		threat_covid safety_perception 										/// controls: perception COVID-19
		school_system_3  medical_system_2 medical_system_1  religion 		/// controls: attitudes towards esoteric and medical belief systems & religion
		left_right GAL_TAN sonntagsfrage spaltungsindex						/// controls: political attitudes  
		solidarity_index big_5_e big_5_n big_5_con big_5_com big_5_o  		/// controls: psychosocial 
		bundesland gender age_group_r edu_3 hh_inc  health_01   			/// controls: sociodemographic     
		informed_vaccination private_health_insurance vacc_reason_* 	if filter_1 == 1
		
	esttab using "tables\Table_S08a.rtf", cells("count mean(fmt(2)) sd(fmt(2)) ") noobs label replace 

	clear all
			
	* Read data 
	use "dump\PPD_July_2022_labeled.dta", replace			
	
	asdoc 	tab framing , replace save(dump\table_S08b_not_unicode.rtf) title(percentages original data) 
	asdoc 	tab framing if filter_1 == 1, append title(percentages listwise deletion) 
	
	asdoc 	tab religion , append  title(percentages original data) 
	asdoc 	tab religion if filter_1 == 1, append title(percentages listwise deletion) 
	
	asdoc 	tab vacc_status , append  title(percentages original data) 
	asdoc 	tab vacc_status if filter_1 == 1, append title(percentages listwise deletion)	
	
	asdoc 	tab sonntagsfrage , append  title(percentages original data) 
	asdoc 	tab sonntagsfrage if filter_1 == 1, append title(percentages listwise deletion)	
	
	asdoc 	tab bundesland , append  title(percentages original data) 
	asdoc 	tab bundesland if filter_1 == 1, append title(percentages listwise deletion) 
	
	asdoc 	tab gender , append  title(percentages original data) 
	asdoc 	tab gender if filter_1 == 1, append title(percentages listwise deletion)	
	
	asdoc 	tab age_group_r , append  title(percentages original data) 
	asdoc 	tab age_group_r if filter_1 == 1, append title(percentages listwise deletion)	
	
	asdoc 	tab edu_3 , append  title(percentages original data) 
	asdoc 	tab edu_3 if filter_1 == 1, append title(percentages listwise deletion) 
	
	asdoc 	tab hh_inc , append  title(percentages original data) 
	asdoc 	tab hh_inc if filter_1 == 1, append title(percentages listwise deletion)		

	asdoc 	tab private , append  title(percentages original data) 
	asdoc 	tab private if filter_1 == 1, append title(percentages listwise deletion)	
	
	unicode convertfile dump\table_S08b_not_unicode.rtf  tables\table_S08b.rtf, dstencoding(Windows-1252) replace
	
	
